Convergence Acceleration via Combined Nonlinear-Condensation Transformations 
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A method of numerically evaluating slowly convergent 
monotone series is described. First, we apply a condensa- 
tion transformation due to Van Wijngaarden to the original 
series. This transforms the original monotone series into an 
alternating series. In the second step, the convergence of the 
transformed series is accelerated with the help of suitable non- 
linear sequence transformations that are known to be partic- 
ularly powerful for alternating series. Some theoretical as- 
pects of our approach are discussed. The efficiency, numerical 
stability, and wide applicability of the combined nonlinear- 
condensation transformation is illustrated by a number of ex- 
amples. We discuss the evaluation of special functions close 
to or on the boundary of the circle of convergence, even in 
the vicinity of singularities. We also consider a series of prod- 
ucts of spherical Bessel functions, which serves as a model for 
partial wave expansions occurring in quantum electrodynamic 
bound state calculations. 
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I. INTRODUCTION 

Divergent and slowly convergent series occur abun- 
dantly in the mathematical and physical sciences. Ac- 
cordingly, there is an extensive literature on numerical 
techniques which convert a divergent or slowly conver- 
gent series into a new series with hopefully better nu- 
merical properties. An overview of the existing sequence 
transformations as well as many references can be found 
in books by Wimp |lj and Brezinski and Redivo Zaglia 
Q . The historical development of these techniques up to 
1945 is described in a monograph by Brezinski |^], and 
the more recent developments are discussed in an article 
by Brezinski ||. 

A very important class of sequences {s n }^_ is char- 
acterized by the asymptotic condition 



lim 

n — >oo s 



Sn+1 



p. 



(1.1) 



which closely resembles the well known ratio test for in- 
finite series. Here, s — Soo is the limit of this sequence 
as n — > oo. A convergent sequence satisfying condition 



(1.1) with \p\ < 1 is called linearly convergent, and it is 
called logarithmically convergent if p = 1. 



If the elements of the sequence {s n }^L m Eq. (1-1) 
are the partial sums s n = X)fc=o ak °^ an infinite series, 
and if p satisfies either p=lorO< J o<l (these are 
the only cases which will be considered in this article), 
then there exists an integer N > such that all terms 
ak with k > N have the same sign. Hence, series of 
this kind can be split up into a finite sum containing the 
leading terms with the irregular signs, and a monotone 
series whose terms all have the same sign. The subject of 
this article is the efficient and reliable evaluation of series 
that exhibit these properties. 

The partial sums of a nonterminating Gaussian hyper- 
geometric series 
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where (a) m = T(a + m)/r(a) = a(a + 1) . . . (a + m — 1) 
is a Pochhammer symbol, or its generalization 
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which converge for \z\ < 1 and diverge for \z\ > 1, are 
typical examples of linearly convergent sequences with 
p = z. 

The partial sums of the Dirichlet series 



C(z) 
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m=0 



(TO + 1)" 



(1.4) 



for the Riemann zeta function converge logarithmically 
if Re(z) > 1. This follows from the following asymptotic 
estimate (n — > oo) of the truncation error (p. 21 of Ref. 
0): 
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The Dirichlet series (1.4) is notorious for its extremely 
slow convergence if Re(z) is only slightly larger than one. 
In this case, the series can only be used for the compu- 
tation of £(z) if it is combined with suitable convergence 
acceleration methods like the Euler-Maclaurin summa- 
tion (see for instance Chapter 8 of Ref. ||, p. 379 of Ref. 
§, or Chapter 6 of Ref. §). 

The acceleration of linear convergence is comparatively 
simple, bot h theoretically and practically, as long as p 
in Eq. (1.1) is not too close to one. With the help of 
Germain-Bonne's formal theory of convergence accelera- 
tion |J and its extension (Section 12 of Ref. 0), it can be 
decided whether a sequence transformation is capable of 
accelerating linear convergence or not. Moreover, many 
sequence transformations are known which are capable 
of accelerating linear convergence effectively. Examples 
arc the A 2 process, which is usually attributed to Aitken 
although it is in fact much older (p. 90 of Ref. 
Wynn's epsilon algorithm jll| , which produces Pade ap- 
proximants if the input data are the partial sums of a 
formal power series, or Levin's sequence transformation 
p2[ and generalizations (Sections 7-9 of Ref. Q), which 
require as input data not only the elements of the se- 
quence to be transformed but also explicit estimates for 
the corresponding truncation errors. 

The acceleration of logarithmic convergence is much 
more difficult than the acceleration of linear convergence. 



Delahaye and Germain-Bonne [1J] showed that no se- 
quence transformation can exist which is able to accel- 
erate the convergence of all logarithmically convergent 
sequences. Consequently, in the case of logarithmic con- 
vergence the success of a convergence acceleration pro- 
cess cannot be guaranteed unless additional information 
is available. Also, an analogue of Germain-Bonne's for- 
mal theory of linear convergence acceleration || cannot 
exist. 

In spite of these complications, many sequence trans- 
formations are known which work reasonably well at least 
for suitably restricted subsets of the class of logarithmi- 
cally convergent sequences. Examples are Richardson ex- 
trapolation |1J], Wynn's rho algorithm jl5| and its itera- 
tion (Section 6 of Ref. 0]) as well as Osada's modification 
of the rho algorithm |]16| , Brezinski's theta algorithm |T^ ] 
and its iteration (Section 10 of Ref. j|), Levin's u and 
v transformations p2] and related transformations (Sec- 
tions 7 - 9 of Ref. ||] ) , and the modification of the A 2 
process by Bj0rstad, Dahlquist, and Grosse Q. How- 
ever, there is a considerable amount of theoretical and 
empirical evidence that sequence transformations are in 
general less effective in the case of logarithmic conver- 
gence than in the case of linear convergence. 

Numerical stability is a very important issue. A se- 
quence transformation can only accelerate convergence 
if it succeeds in extracting some additional information 
about the index-dependence of the truncation errors from 
a finite set s n , s n+ \, . . ., s n+ k of input data. Nor- 
mally, this is done by forming arithmetic expressions in- 
volving higher weighted differences. However, forming 
higher weighted differences is a potentially unstable pro- 
cess which can easily lead to a serious loss of significant 
digits or even to completely nonsensical results. 

If the input data are the partial sums of a strictly al- 
ternating series, the formation of higher weighted differ- 
ences is normally a remarkably stable process. Hence, a 
serious loss of significant digits is not to be expected if 
the partial sums of a strictly alternating series are used 
as input data in a convergence acceleration or summa- 
tion process. If, however, the input data are the partial 
sums of a monotone series, numerical instabilities due to 
cancellation are quite likely, in particular if convergence 
is very slow. Thus, if the sequence to be transf orm ed ei- 
ther converges linearly with a value of p in Eq. (1.1) that 
is only slightly smaller than one, or if it converges loga- 
rithmically (p — 1), numerical instabilities are a serious 
problem and at least some loss of significant digits is to 
be expected. 

Generally, the sequence transformations mentioned 
above are not able to determine the limit of a logarithmi- 
cally convergent series, whose terms ultimately all have 
the same sign, with an accuracy close to machine accu- 
racy. This restricts the practical usefulness of sequence 
transformations severely, e.g., in FORTRAN calculations 
with a fixed precision. 

In this paper, we show that these stability problems 
can often be overcome by transforming slowly convergent 
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monotone series not by straightforward application of a 
single sequence transformation but by a combination of 
two different transformations. 

In the first step, a monotone series is transformed into 
a strictly alternating series with the help of a condensa- 
tion transformation. This transformation was first men- 
tioned on p. 126 of Rcf. (l9) and only later published 
by Van Wijngaarden J2(J. Later, the Van Wijngaarden 
transformation was studied by Daniel and recently 
it was rederived by Pelzl and King |22]| , who used it for 
the high-precision calculation of atomic three-electron in- 
teraction integrals of explicitly correlated wave functions. 

In the second step, the convergence of the resulting al- 
ternating series is accelerated by suitable nonlinear se- 
quence transformations f^,[^| which are known to be 
very powerful in the case of alternating series. Since 
the transformation of alternating series is a remarkably 
stable process, the limits of even extremely slowly con- 
vergent monotone series can be determined with an ac- 
curacy close to machine accuracy. Conceptually, but not 
technically our approach resembles that of Brezinski, De- 
lahaye, and Germain-Bonne p3j who proposed to extract 
a linearly convergent subsequence from a logarithmically 
convergent input sequence by a selection process. 

In this article, we will call our approach, which con- 
sists of the Van Wijngaarden condensation transforma- 
tion and the subsequent nonlinear sequence transforma- 
tion, the combined nonlinear- condensation transforma- 
tion (CNCT). 

The CNCT is not restricted to logarithmically conver- 
gent series. It can also be used in the case of a lin early 
convergent monotone series with a value of p in Eq. (1.1) 
close to one. Typically, this corresponds to a power series 
whose coefficients ultimately all have the same sign and 
whose argument is positive and close to the boundary of 
the circle of convergence. We will present some examples 
which show that the CNCT works even if the argument 
of the power series is very close to a singularity. 

Our approach requires the evaluation of terms of the 
original monotone series with high indices. Consequently, 
our two-step approach is computationally more demand- 
ing than the application of a single sequence transfor- 
mation, and it cannot be applied if only a few terms 
of a slowly convergent series are available. In spite of 
these restrictions, we believe that the CNCT is very use- 
ful at least for certain problems since it is able to produce 
highly accurate results that can only be accomplished 
otherwise with a considerably greater numerical effort. 

Special functions are defined and, in many cases, also 
evaluated via series expansions. The evaluation of special 
functions is an old problem of numerical mathematics 
with a very extensive literature (compare for example 
the books by Luke 24 25fl, Van Der Laan and Temme 
fjrjf , and Zhang and Jin p7|). Nevertheless, there is still 
a considerable amount of research going on and many 
new algorithms for the computation of special functions 
have been developed recently (compare for example the 
papers by Lozier and Olver |28] and Lozier f29] and the 



long lists of references therein). We believe that sequence 
transformations in general and the CNCT in particular 
are useful tools for the evaluation of special functions. Of 
course, this is also true for problems in theoretical and 
computational physics. 

This paper is organized as follows: In Section [nj, we de- 
scribe the Van Wijngaarden transformation. In Section 



[II 



we discuss the nonlinear sequence transformations 
which we use for the acceleration of the resulting alter- 
nating series. In Section [V, we apply the CNCT to the 
Riemann zeta function. In Section we consider the 
evaluation of the Lerch transcendent and related func- 
tions with arguments on or close to the boundary of the 
circle of convergence. In Section yj\ we examine the 
evaluation of the generalized hypergeometric series p +\F p 
(p > 2) with unit argument or with an argu ment z which 
is only slightly smaller than one. In Section VII , we dis- 
cuss the evaluation of an infinite series involving Bessel 
and Hankel functions. In Appendix |a|, we discuss the ef- 
ficiency of sequence transformations in the case of mono- 
tone series or strictly alternating series in more detail. 
Finally, in Appendix [B] we discuss exactness properties 
of the sequence transformations which we use in the sec- 
ond step for the acceleration of the convergence of the 
resulting alternating series. 



The example of Section |VI]| serves as a model prob- 
lem for slowly convergent series, which occur in quantum 
electrodynamic bound state calculations and which were 
treated successfully with the methods discussed in this 
paper [jjo). Therefore, we expect the CNCT to be a gen- 
eral computational tool for the evaluation of slowly con- 
vergent sums over intermediate angular momenta which 
arise from the decomposition of relativistic propagators 
in QED bound state calculations into partial waves. 

All calculations were done in Mathematical with a rel- 
ative accuracy of 16 decimal digits pi]] . In this way, we 
simulate the usual DOUBLE PRECISION accuracy in 
FORTRAN. 



II. THE VAN WIJNGAARDEN 
TRANSFORMATION 



Let us assume that the partial sums 

n 

cr„ = ^ a ( k ) 



(2.1) 



k=0 



1 Certain commercial equipment, instruments, or materials 
are identified in this paper to foster understanding. Such 
identification does not imply recommendation or endorsement 
by the National Institute of Standards and Technology, nor 
does it imply that the materials or equipment identified are 
necessarily the best available for the purpose. 
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of an infinite series converge either linearly or logarithmi- 
cally to some limit a = as n — ► oo. We also assume 
that all terms a(k) have the same sign, i.e., the series 
^2kLo a (k) i s a monotone series. 

Following Van Wijngaarden p(J, we transform the 
original series into an alternating series 53fc=o( — -0 3 'Ai 
according to 



^o(fe)=x;(-i) i A i , 



fe=0 



fc=0 



,0) 



= 2 fc a(2 fc (. ? + !)-!) 



(2.2) 

(2.3) 
(2.4) 



Obviously, the terms Aj defined in Eq. ( [2.3j ) all have 
the same sign if the terms a(k) of the original series all 
have the same sign. In the sequel, the quantities Aj 
will be referred to as the condensed series, and the series 
^2j =0 (— 1) J Aj will be referred to as the transformed al- 
ternating series, or alternatively as the Van Wijngaarden 
transformed series. 

We call the Van Wijngaarden transformation a con- 
densation transformation because its close connection to 
Cauchy's condensation theorem (p. 28 of Ref. [B2j or p. 
121 of Ref. P3|| ). Given a monotone series Y^kLo 11 ^) 
with terms that satisfy \a(k + 1)| < |a(fc)|, Cauchy's con- 
densation theorem states that X)fc=o a (^) converges if and 
only if the first condensed series Ao defined in Eq. (2.5) 
converges. 

In Eq. ( |2.3| ), the indices of the terms of the original 
series are chosen in such a way that sa mpl ing at very high 
indices takes place (according to Eq. (2.4), the indices of 
the terms of the original series grow exponentially). In 
this way, we obtain information about the behaviour of 
the terms of the original series at high indices. 

Moreover, if the terms of the original series behave 
asymptotically (n — > oo) either like a(n) ~ n~ 1 ~ £ with 
e > or like a(n) ~ n"r n with < r < 1 and (3 real, then 
the terms of the original series become negligibly small 
after a few evaluations. Specifically, the series ( [2.3] ) for 
the terms Aj converges linearly in t hese cases. 

When summing over k in Eq. ( |2.3| ), we found that it is 
normally sufficient to terminate this sum when the last 
term is smaller than the desired accuracy. Typically, 20 
to 30 terms are needed for a relative accuracy of 10~ 14 
in the final result for the condensed sum. 

It shoul d be possible to accelerate the convergence of 
the series (2.3) for Aj. Since, however, this is a monotone 
series, it is not clear whether and how many digits would 
be lost in the convergence acceleration process. Conse- 
quently, we prefer to perform a safe and straightforward 
evaluation of the condensed sums Aj and add up the 
terms of the series until convergence is reached, although 
this is most likely not the most efficient approach. 

The transformation from a monotone series to a 



essentially a reordering of the terms a(k) of the original 
series. This is seen as follows. We first define the partial 
sums 



E 

j=0 



(2.5) 



of the Van Wijngaarden transformed original series. It 
can be shown easily that S n with n > reproduces the 
partial sum a n , Eq. ( [2.1[ ), which contains the first n + 1 
terms of the original series. To illustrate this procedure, 
we present Table |. Formal proofs of the correctness of 
this rearrangement can be found in Ref. J2jj ] or in the 
Appendix of Ref. @. 

Insert Table || here 

Daniel |nj was able to formulate some mild conditions 
which guarantee that the limits a and S of the partial 



and are equal: 
lim 



sums (2.1) and (2.5), respectively, simultaneously exist 



= a = lim S„ — S . 



(2.6) 



For example, in the Corollary on p. 92 of Ref. |2l| it was 
shown that if a strictly decreasing sequence {Mk}"j!° =0 of 
positive bounds exists which satisfy \a(k)\ < Mk for all 
k > and if X^fclo ^ k ^ 00 holds, then the original 
monotone series and the Van Wijng aard en transformed 
series on the right-hand side of Eq. (2.2) both converge 
to the same limit (i.e., a = S holds). This useful criterion 
is fulfilled by all series considered in this paper. 



III. NONLINEAR SEQUENCE 
TRANSFORMATIONS 



The series (2.3) for the terms of the Van Wijngaarden 
transformed series can be rewritten as follows: 



Aj = a(j) + 2a(2j + 1) + 4a(4j + 3) + 



(3.1) 



Since the terms a(fc) of the original series have by as- 
sumption the same sign, we immediately observe 



I A, I > la,- 



(3.2) 



strictly alternating series according to Eqs. (2J2) - (2A) is 



Con seq uent ly, the Van Wijngaarden transformation, Eqs. 
( |2.2D - ( |2.4| ) , does not lead to an alternating series whose 
terms decay more rapidly in magnitude than the terms 
of the original monotone series. Thus, an acceleration 
of convergence can only be achieved if the partial sums 
( |2.5D of the Van Wijngaarden transformed series are used 
as input data in a convergence acceleration process. 

Since the Van Wijnga ard en transformed series on the 
right-hand side of Eq. ( |2.2| ) is alternating if the terms 
of the original series all have the same sign, it is rec- 
ommended to choose a suitable convergence acceleration 
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method which is particularly powerful in the case of al- 
ternating series. A judicious choice of the convergence 
accelerator is of utmost importance since it ultimately 
decides whether our approach is numerically useful or 
not. 



Daniel 21 used the Euler transformation (Eq. (3.6.27) 
on p. 16 of Ref. §1)): 



CO CO 

£(-D fc «* = £ 



k=0 



k=0 



(-1)* 
2 fc+l 



A k u . 



(3.3) 



Here, A is the (forward) difference operator defined by 
A/(») = /(n+l)-/(n), and 



A k u = (-l) fe J2 



m=0 



(3.4) 



The Euler transformation, which was published in its 
original version in 1755 on p. 281 of Ref. J35|, is a se- 
ries transformation which is specially designed for alter- 
nating series. It is treated in many books on numerical 
mathematics. Nevertheless, the Euler transformation is 
not a particularly efficient accelerator for the Van Wijn- 
gaarden transformed series considered in this article. In 
Table H we show some explicit results obtained by the 
Euler transformation. We also applied the Euler trans- 
formation to all other Van Wijngaarden transformed se- 
ries presented in this paper, and we consistently observed 
that it is clearly less powerful than the transformations 
we used. 

Much better results can be expected from the more 
modern nonlinear sequence transformations which trans- 
form a sequence {sn}^L ' wnose elements may be the 
partial sums of an infinite series, into a new sequence 
{ s 'n}ri=o w ith better numerical properties [[l]||||. 

The basic assumption of a sequence transformation is 
that the elements of the sequence {s„}™ =0 to be trans- 
formed can be partitioned into a limit s and a remainder 
r n according to 

s„ = s + r n (3.5) 

for all n > 0. Only in the case of some more or less triv- 
ial model cases will a sequence transformation be able 
to determine the limit s of the sequence {s n }^L exactly 
after a finite number of steps. Hence, the elements of the 
transformed sequence {s^}^L can be partitioned into 
the same limit s and a transformed remainder r' n accord- 
ing to 



= s + r 



(3.6) 



for all n > 0. In general, the transformed remainders r' n 
are nonzero for all finite values of n. However, conver- 
gence is obviously increased if the transformed remain- 
ders {^}^Lo vanish more rapidly than the original re- 
mainders {r n }^Lo- 



fn+l ,. Sn+l 

lim = lim 



= 0. 



n— too TV 



n— too Sr, — S 



(3.7) 



The best known nonlinear sequence transformation is 
probably Wynn's epsilon algorithm |ll]] which produces 
Pade approximants if the input data are the partial sums 
of a formal power series. The epsilon algorithm is also 
known to work well in the case of alternating series. Con- 
sequently, it is an obvious idea to use the epsilon algo- 
rithm for the acceleration of the alternating series which 
we obtain via the Van Wijngaarden transformation. 

However, much better results can be obtained by ap- 
plying sequence transformations which use explicit re- 
mainder estimates as input data in addition to the par- 
tial sums of the series to be transformed. Consequently, 
in this article we use exclusively the following two se- 
quence transformations by Levin p"2fl and Weniger ||, 
respectively: 



£fe (A s n ,u n ) 



£(-d 

3=0 



_ A k {{n + p) k - 1 s n /co n } 
A k {(n + P) k ~ 1 / ljJn } 

k\ (p + n + j)*- 1 s n+J 



EM) 

3=0 



fe-1 



(P + n + j) 

fc-i 



(3.8) 



jj {d + n + kf- 1 u n+j 



Mn) fR * _ A fc {(n + /3) fc -i s n /uj n } 

°k \P,S n ,Wn) — — XT77 ~lr\ 7 T — 

A*- {{n + f3)k-i/Lu n } 

k 
3=0 



j J (p + n + k)k-i u-n+j 



k 

E 

3=0 



(p + n + j) k -i 



(3.9) 



j J (p + n + k)k-i u-n+j 



Here, {s n }^Lo * s the sequence to be transformed, and 
{ w n}^o ls a sequence of truncation error estimates. The 
shift parameter B has to be positive in order to permit 
n = in Eqs. ( |3.8D and ( |3 . 9| ) . The most obvious choice is 
= 1. This has been the choice in virtually all previous 
applications of these sequence transformations, and we 
will use P = 1 unless explicitly stated. However, we 
show in Appendix |B] that in some cases it may be very 
advantageous to choose other values for p. 



T he n umerator and denominator sums in Eqs. (3.8) 
and (3.9) can also be computed with the help of the three- 
term recursions (Eq. (7.2-8) of Ref. M) 



(P + n + k)(P + n + k) k 



fc-i 



and (Eq. 



(P + n + k 
.3-7) of Ref. [I 



1) A 



4%) 



(3.10) 



sl%m = si n+1 \p) 

(P + n + k)(P + n + k-l) 



(P + n + 2k)(P + n + 2k-l) k Kf ' 



(3.11) 
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The initial values 



4%) 



4 n) (/3) = S< n) (/3) = s n /u n and 
Sq U \p) = l/u>n produce the numerator and 
denominator sums, respectively, of C^jf 1 s n , u) n ) and 



The performance of the sequence transformations 
i™ ) (/3,s„,w„) and ' 



depends crucially on 



the remainder estimates {u„}J° =(] . Under exceptionally 
favourable circumstances it may be possible to construct 
explicit approximations to the truncation errors of the 
input sequence {sn}^o- In most cases of practical inter- 
est this is not possible, since only the numerical values of 
a finite number of sequence elements is available. Hence, 
in practice one has to determine the remainder estimates 
from these numerical values. 

On the basis of purely heuristic and asymptotic argu- 
ments, Levin fl2]| proposed some simple remainder esti- 
mates which - when used in Eq. ( |3.8| ) - lead to Levin's 
it, t, and v transformations. These remainder estimates 
can also be used in Eq. (Section 8.4 of Ref. J|). 

However, in the case of a strictly alternating series the 
best simple truncation error estimate is the first term 
neglected in the partial sum (p. 259 of Ref. |33|]). This is 
also the best simple estimate for the truncation error of a 
strictly alternating nontcrminating hypergcometric series 
2-Fo (a, /3; —z) with a,/3,z > 0, which diverges strongly for 
every z^O (Theorem 5.12-5 of Ref. Pq| ). Consequently, 
in the case of convergent or divergent alternating series 
it is a natural idea to use the first term neglected in the 
partial sum as the remainder estimate, as proposed by 
Smith and Ford §7|. 

In this article, we always accelerate the convergence of 
the par tial sums S„ of the Van Wijngaarden transformed 
series (2.3) which is strictly alternating if the original 
series is monoton e. T hus, w e us e exclusively the sequence 
transformations (3.8) and (3.9) in combination with the 
remainder estimate proposed by Smith and Ford p%: 



= AS„ = (-1) 



n+l 



(3.12) 



This yields the following sequence transformations (Eqs. 
(7.3-9) and (8.4-4) of Ref. [§): 



4 n) (/3,S n ) = 4" , (/3,S„,AS„) 



_ r (n) 



4 n) (/3,S„,(-l)" +i A„ +1 ) 



n+l) , 



4 n) (AS, 



4 n) (/3,S„,AS„) 



n+l , 



(3.13) 
(3.14) 



Unless explicitly stated, we use (3—1. In the applica- 
tions described in the following, it will be clear from the 
context which monotone series was transformed accord- 
ing to Van Wijngaarden such as to produce the input 
data for the nonlinear sequence transformations which 
are the partial sums S„ of the alternating series (2.3). 
Alternative remainder estimates for the sequence 
and ( |3.9| ) were discussed in Sec- 
of ReT § or in Refs. ||,||. 



transformations (3 
tions 7 and 



j, From a purely theoretical point of view, the sequence 



transformations d k and S k l> as well as their parent 



r(n) 



transformations C k n ^ and S k ' l> have the disadvantage that 
no general convergence proof is known. Only for some 
special model problems could rigorous convergence proofs 
be obtained (Refs. @-£|| or Sections 13 and 14 of Ref. 
P|). However, there is overwhelming empirical evidence 



:(") 



that and S [ k n> work very well in the case of conver- 
gent or divergent alternating series for instance as they 
occur in special function theory ||||[l^,[37 43- 48|] or in 
quantum mechanical perturbation theory H IM , [§! 49 - 55|] . 

Pelzl and King [£2| oru y used the Levin transforma- 
tion d^ S n ) for the acceleration of the convergence of 
Van Wijngaarden transformed alternating series. How- 
ever, we shall show that the closely related transforma- 
tion (/?, S„) frequently gives better results. 



I 



Application of the sequence transformations (3.13) and 
14) to the partial sums S n of the infinite series (2.3) 
produces doubly indexed sets of transforms, which de- 
pend on the starting index n of the input data and the 
transformation order k: 



(3.15) 



Here, T ( k n) stands for either d{ n) (/3, S„) or 8 ( k n) (/3, S„) . 

The transforms can be displayed in a rectangular 
scheme called the table of the transformation. In prin- 
ciple, there is an unlimited variety of different paths, on 
which we can move through the table in a convergence 
acceleration or summation process. In this article, we 
always proceed on such a path that for a given set of 
input data, we use the transforms with the highest possi- 
ble transformation order as approximations to the limit 
of the series to be transformed. Thus, given the input 
data {So, Si, . . . , S m , S m +i} we always use the sequence 
| j,(o) ^ ji(o) ^ ji(fp | Q £ t rans f orms ag approximations to 

the limit (compare for instance Eqs. (7-5.4) and (7-5.5) 

of Ref. H ) , where Tm^ should provide the best approxi- 
mation since it has the highest transformation order. 

In Appendix we explain why the transformations 
Cf? (f3, s n ,u> n ) and Sj^' (ft, s n ,u) n ) are much more effec- 
tive in the case of an alternating series than in the case 
of a monotone series. 



IV. THE RIEMANN ZETA FUNCTION 

The Riemann zeta function is discussed in most books 
on special functions. Particularly detailed treatments can 
be found in monographs by Edwards Q , Titchmarsh [|56f , 
Ivic ^7j, and Patterson ]58| ]. Many applications of the 
zeta function in theoretical physics are described in books 
by Elizalde, Odintsov, Romeo, Bytsenko, and Zerbini p9| ] 
and by Elizalde @. 

In this section, we discuss how the CNCT can be used 
for the evaluation of the Riemann zeta function, start- 



G 



ing from the logarithmically convergent Dirichlet series 
(1.4). We do not claim that our approach is necessarily 
more powerful than other, more specialized techniques 
for the evaluation of the Riemann zeta function. How- 
ever, because of its simplicity, the Dirichlet series ( [Oh 
is particularly well suited for an illustration of our ap- 
proach. 

In the case of the Riemann zeta function, the terms 



bj? of the series (2.3) for Aj are given by 



.0) 



(2- 



U + 1) 2 



(4.1) 



Thus, the series (2J3) is a power series in 2 1 ~ z which con- 
verges linearly. Moreover, it is essentially the geometric 
series in the variable 2 1_z that can be expressed in closed 
form according to 



A, 



1 



1 



0' + l Y 



(4.2) 



Inserting th is in to the infinite series on the right-hand 
side of Eq. (2.2) yields the following transformed alter- 
nating series for the Riemann zeta function: 



C(*) = 



1 



1 



E 



(4.3) 



for this argument, but its convergence is so slow that it 
is computationally useless. It follows from the trunca- 
tion error estimate ( [L.5| ) that on the order of 10 600 terms 
would be needed to achieve the modest relative accuracy 
of 10~ 6 if the Dirichlet series were summed by adding up 
the terms. A much more efficient evaluation of £(z) near 
Re(z) = 1 is based on the Euler-Maclaurin formula (Eqs. 
(2.01) and (2.02) on p. 285 of Ref. @) 

E /(") = / H*)** + oL/W + ZW] 

n=N jN 1 



Ra 



i 



/W-i)(A0-/ (2j '- 1) Wl + R g , (4.4a) 



N 



B 2q (x-lx})f^(x)dx. (4.4b) 



Here, [x] is the integral part of x, i.e., the largest integer 
m satisfying m < x, Bk(x) is a Bernoulli polynomial 
defined by the generating function (Eq. (1.06) on p. 281 
of Ref. K) 



texp(a;i) 
exp(t) — 1 



E B ^ 3 



\t\ < 2tt. 



(4.5) 



and 



An alternative proof of the validity of this series expan- 
sion can be found in Section 2.2 of Titchmarsh's book 



The terms of the alternating series (4.3) decay as slowly 
in magnitude as the ter ms of the logarithmically con- 
vergent Dirichlet series (1.4). Nevertheless, the series 
(4.3) offers some substantial computational advantages 
because it is alternating. For example, it converges for 
Re(z) > 0, whereas the Dirichlet series converges onl y fo r 
Re(z) > 1. For Re(z) < both the Dirichlet series flO| ) 
and the alternating series (13) diverge. However, the al- 
ternating series can provide an analytic continuation of 
the Riemann zeta function for Re(z) < if it is combined 
with a suitable summation process. For this p urpos e, we 
can use the same nonlinear transformations (3.13) and 
(| 3.14 ) as in the case of a convergent Van Wijngaarden 



transformed alternating series. 

This observation extends the applicability of the 
CNCT, which was originally designed for convergent 
monotone series only. Of course, this extension to diver- 
gent series is only possible if the monotone series (2.2) 



for Aj can be summed. Normally, the summation of a 
monotone divergent series is very difficult. However, in 
the case of the Riemann zeta function this is trivial since 
the series (2.3) is - as remarked above - essentially the 
geometric series in the variable 2 1-2 which can be ex- 



pressed in closed form according to Eq. (4.2), no matter 
whether (2 1 -*! < 1 or \2 X - Z \ > 1 holds. 

Our first numerical example is the zet a fun ction with 
argument z = 1.01. The Dirichlet series ( |l.4| ) converges 



Br, 



B m (0) 



(4.6) 



is a Bernoulli number (p. 281 of Ref. j5|). 

In the next step, we rewrite the Dirichlet series (|l.4| ) 
as follows, 



C(z) 



1 



N-l 

E ( m + 1)2 

m=0 v ; 



E 



-N 



(m + l) 2 



(4.7) 



and apply the Euler-Maclaurin formula (4.4) to the in- 
finite series on the right-hand side. With N = 100, we 
easily obtain the value of the zeta function accurate to 
more than 15 significant decimal digits. The value of 
£(1.01) to 15 decimal places is 



10~ J C(1-01) = 0.100 577 943 338 497. 



(4.8) 



We present this result (as well as all other subsequent 
numerical results) normalized to a number in the interval 
(0, 1) by a multiplicative power of 10. 

In Table O, we apply the Euler transform ation (|3.3| ) 
and t he nonlinear sequence transformations (3.13) and 
(|3 .14 ) to the partial sums of the alternating series (|4.3| ) 
with z — 1.01. We list as a function of n the partial sums 



E 



(-i) fe 



1-2 1 - 2 ^ (fc + 1) 
fc=0 v ; 



(4.9) 



of the alternating series (4.3), the partial sums 
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E, 



k (-1) 
1 x ^ 1 i. -\ \m 

I _ 2 1 - 2 2' 1 '+ 1 



fc=0 m=0 



(m + 1) 2 



(4.10) 



of the Euler transformed series ( p.3[ ) , and the nonlinear 
sequence transformations (l, So) and (l, So) ■ 

Insert Table |T| here 

The results in Table |j| show that the Euler transfor- 
mation is in the case of the alternating series (L3) much 
less effective than the nonlinear sequence transformations 

dn (l, S ) and 5 ^ ( l, S ) . When we applied Wynn's ep- 
silon algorithm Jul to the alternating series (4.3), we 
also obtained clearly inferior results. This is in agree- 
ment with the observations of Pelzl and King |22| . 

Thus, in the follow ing w e will only c onsider the nonlin- 
ear transformations (^13|) and ( [3T4l) which use explicit 
remainder estimates. 

In Table ||, we present apparently redundant data 
since our transformation results in the last two columns 
do not change for n > 12. However, we include these re- 
sults here as well as below in order to demonstrate that 
our transformation remains stable even if we increase the 
transformation order. 

If the argument z is zero or a negative integer, the 
Riemann zeta function satisfies (p. 807 of Ref. Q): 



C(o) 

C(l - 2m) 



1 

2 ' 

-E>2m 

2m 



0. 



C(-2m) 

m = 1,2 



(4.11a) 
(4.11b) 



Here, i?2m is a Bernoulli number defined in Eq. ( |4.6| ). 

Our second numerical example is the zeta function 
with argument z = — 1 which, because B2 = 1/6, sat- 
isfies 



C(-i) 



-1/12 



-0.0833 333. 



(4.12) 



As remarked above, the alternating series (4.3) diverges 
for z = — 1. However, our results in Table III show that 



this series can be summe d effe ctivel y by the nonlinear 
sequence transformations (3.13) and (3.14). 



Insert Table |TI| here 



The results in Table [II indicate that the sequence 
transfo rma tion <J„ (l, So) is exact for the alternating 
series ( |4.3[ ) with z = —1. The exactness of the se- 

and 



quence transformations C^\/3 



,u n ), Eq. (3 



uj n ), Eq. (f^Sl), for the infinite series (|4.3|) and 



for related series is discussed in Appendix [Bj. 

There is considerable research going on in connec- 
tion with the Riemann-Siegel conjecture |n] pi| . In this 



context, it is necessary to evaluate the Riemann zeta 
function on the so-called critical line z = 1/2 + t i 
(—00 < t < 00). This can also be a ccom plishe d by ap- 
plying the sequence transformations ( 3. IS ) an d ( 3.14 ) to 
the partial sums of the alternating series (4.2). In Table 
IV, where we treat the case z = 1/2+13. 7 i, the Levin 
transformation outperforms the Weniger transformation. 
For limitations of space, we do not present the partial 
sums of the transformed alternating series in Table IV. 



Insert Table [TV] here 

V. THE LERCH TRANSCENDENT AND 
RELATED FUNCTIONS 

The Lerch transcendent (p. 32 of Ref. p5[) 

00 n 

®(z,s,a) = 1 1 \ z \< 1 > 
* — ' (a + n) s 

contains many special functions as special cases, for ex- 
ample the Riemann zeta function 



C(s) = 9(1,8,1) 
or the Jonquiere function (p. 33 of Ref. |6 

n+l 



(n + l) s 



= z$(z,s,l). 



(5.2) 



(5.3) 



In the physics literature, the Jonquiere function is usually 
called generalized logarithm or polylogarithm, and the 
following notation is used: 



Li s (z) = £ 



yk+l 



k=0 



(k + iy 



F(z,s) 



(5.4) 



This is also the notation and terminology used in the 
book by Lewin |36| . 

Lerch transcendents and their special cases, the poly- 
logarithms, are ubiquitous in theoretical physics. They 
play a role in Bose-Einstein condensation fS7|]6"g| l, and 
they are particularly important in quantum field theory. 
For example, they occur in integrands of one-dimensional 
numerical integrations in quantum electrodynamic bound 
state calculations fp9r|7S]. 



The series expansion (5.1) for Q(z, s, a) converges very 
slowly if the argument z is only slightly smaller than 
one. However, the CNCT can be used for an efficient 
and reliable evaluation of the Lerch transcendent and its 
special cases. 

In our first example (Table [v]) , we evaluate 

10~ 2 Lii (0.99999) = -10~ 2 ln(0.00001) 

= 0.115 129 254 649 702. (5.5) 
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The series (5.4) for Li s (z) converges linearly for \z\ < 1, 
but 



VI. THE GENERALIZED HYPERGEOMETRIC 
SERIES 



x z k+1 



(5.6) 



k=0 



has a singularity at z = 1. Thus, in Table ^ we evaluate 
Lii(z) in the immediate vicinity of a singularity. 

Insert Table [V] here 

In contrast, the series expansions 

yk+l 



U 2 (z) = J2 



fc=0 



and 



Lia(*) = £ 



fc=0 



(fc+1) 2 



,k+l 



(fc+1) 3 



(5.7) 



(5.8) 



converge logarithmically for z = 1. In Tables [v| and VII , 
we consider 

10 _1 Li 2 (0.99999) = 0.164 480 893 699 293... (5.9) 
and 

lO" 1 Li 3 (0.99999) = 0.120 204 045 438 733 , (5.10) 
respectively. As in most previous examples, the Weniger 



transform <5„ (l,So) produces in Tables VI and VII 
somewhat better results than the Levin transform 



Insert Tables [VT] and |V11| here 



In our last example of this section, we evaluate 

10 4 $(0.99999, 2, 10000) = 0.798 585 139 222 548 

(5.11) 

with the help of the CNCT. The prefactor of 10000 is 
introduced so that the final result is of order one. 

Insert Table [VTTT| here 



The Gaussian hypergeometri c fu nction 2 Fi(a, b; c; z), 
defined by the series expansion (1.2) for \z\ < 1, is one of 
the most used and best understood special functions. It 
is discussed in virtually all books on special functions. A 
particularly detailed treatment can be found in the book 
by Slater @. 

The fact that the mathematical properties of the hy- 
pergeometric function 2 F\ are so well u nder stood, greatly 
facilitates its computation. The series ( |l.2| ) does not suf- 
fice for the computation of a nonterminating 2-F1 since it 
converges only if \z\ < 1. By contrast, the hypergeomet- 
ric function 2 F\ is a multivalued function defined in the 
whole complex plane with branch points at z = 1 and 00. 

However, a hypergeometric function 2 Fi{a, b; c; z) can 
be transformed into the sum of two other 2-F1 's with argu- 
ments w = 1 — z, w — 1/z, w — 1/(1 — z), or w = 1 — 1/z, 
respectively. Thus, the argument w of the two resulting 
hypergeometric series can normally be chosen in such a 
way that the two new series in w either converge, if the 
original series in z diverges, or they converge more rapidly 
if the original series converges too slowly to be numeri- 
cally useful. 

For example, if |1 — z\ < 1 and if c — a — b is not a 
positive or negative integer, then we can use the analytic 
continuation formula (Eq. (15.3.6) of Ref. p3) 



2 Fi(a,b; c;z) 
_ T(c)T(c -a-b) 
~ T(c-a)T(c-b) 
T(c)T(a + b- c) 



2-Fi(a, b; a + b — c+ 1; 1 — z) 
^_ z y-a-b 



T(a)T(b) 

x 2-Fi(c — a, c — b; c— a — b + 1; 1 — z) . (6.1) 

Obviously, the two hypergeometric series in 1 — z will 
converge rapidly in the vicinity of z — 1. With the help of 
this or other analytic continuation formulas it is possible 
to compute a hypergeometric function 2-Fi(a ) b; c; z) with 
arbitrary real argument z S (—00, +00) effectively via the 
resulting hyper geo metric series (see p. 127 of Ref. |75) or 
Table I of Ref. Q). 

The situation is much less favourable in the case of the 
generalized hypergeometric series 



i F 2 (a,b, c; d,e;z) 



(a)m(b) m (c) m z m 
(d) m (e) m ml 



(6.2) 



and it gets progressively worse in the case of higher gen- 
eralized hypergeometric series P +\F p with p > 3. The 
analytic continuation formulas of the type of Eq. (6.1) 



are not always known, and those that are known become 
more and more complicated with increasing p [f77| p0[ . 
Thus, the efficient and reliable computation of a gener- 
alized hypergeo metr ic function p +\F p with p > 2 via its 
series expansion HL3) is still a more or less open problem. 
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For example, in Theorem 2 of Ref. [M explicit expres- 
sions for the analytic continuation formulas of a p +\F p 
with p = 2, 3, . . . around z = 1 were constructed. How- 
ever, these expansions are by some orders of magn itude 
more complicated than the analogous formula (6.1) for a 

We show here that the CNCT can be very useful if 
the argument z of a generalized hypergeometric series 
p+\F p is only slightly smaller than one, even if z = 1 is a 
singularity. 

The function 3 F 2 (1, 3/2, 5; 9/8, 47/8; z) has a singular- 
ity at z = 1 because the sum of the real parts of the 
numerator parameters is greater than the sum of the real 
parts of the denominator parameters (see p. 45 of Ref. 
|f74|). Thus, in Table [X we evaluate the function 



1(T 4 3 F 2 (1, 3/2, 5; 9/8, 47/8; 0.99999) 
= 0.238 434 298 763 330 



(6.3) 



in the immediate vicinity of the singularity. With the 
help of the CNCT, it is nevertheless possible to evaluate 
this function without any noticeable loss of significant 
digits. 



Insert Table [1X1 here 



of gamma functions as in Eq. (6.5), but an infinite series 
with complicated terms. 

Otherwise, there is a variety of different simple sum- 
mation theorems for generalized hypergeometric series 
p+\F p with unit argument which are ratios of gamma 
function. However, these summation theorems depend 
upon the value of p, and are usually only valid for cer- 
tain values and/or combinations of the parameters a,\, 
a2, . • • , a p +i and /3±, @2, . . . , p . A list of these sim- 
ple summation theorems can be found in Appendix III 
of Slater's book j74|]. The reconstruction of known sim- 
ple summation theorems for nonterminating generalized 
hypergeometric series and the construction of new ones 
with the help of computer algebra systems is discussed in 
books by Petkovsek, Wilf, and Zeilberger |8l[| and Koepf 

With the help of Watson's summation theorem (p. 245 
Ref. [frl) we obtain 



10" 



5 F 2 (1,3,7;5/2,14;1) = 
0.267 108 047 538 428. 



567567 7T 2 
20971520 



(6.6) 



We show in Table XI that this result can also be obtained 
with the help of the CNCT, starting from the hypergeo- 



metric series (3.2). 



The function 3^2(1, 3, 7; 5/2, 14; z) does not have a sin- 
gularity at z = 1, i.e., the hypergeometric series (|6.2| ) 
converges for z = 1, albeit very slowly. In Table |X| we 
compute the function 



Insert Table IXII here 



10" 1 3 F 2 (1, 3, 7; 5/2, 14; 0.99999) 
= 0.267 102 823 984 762 



(6.4) 



VII. PRODUCTS OF BESSEL AND HANKEL 
FUNCTIONS 



with the help of the CNCT. 

Insert Table |X| here 

There is another major difference between a Gaussian 
hypergeometric series 2 -Pi and a generalized hypergeo- 
metric series p +iF p with p > 2. The series ( |l.2[ ) for a hy- 
pergeometric function 2-^1(0, b; c; z) with unit argument 
z = 1 converges if Re(c) > Re(a+6), and its value is given 
by the Gauss summation theorem (compare for instance 
Section 1.7 of Ref. @): 



2 Fi(a,b;c; 1) = 



r(c)r( c 



b) 



r(c~ a)T(c-b) 



(6.5) 



The series ( [L3j ) for p+1 F p (cti, . . . ,a p+ i; Pi, . . . , /3 P ; z) 
with p > 2 converges at unit argument z = 1 if 
Re(/?i + . . . + ftp) > Rc(cti + . . . + a p +i) (compare p. 
45 of Ref. Q). In Theorem 1 of Ref. [g0|, an analogue 
of the Gauss summation theorem for a generalized hy- 
pergeometric series p +\F p with unit argument was for- 
mulated. However, this expression is not a simple ratio 



In this section, we want to apply the CNCT to series 
whose terms involve more complicated entities and which 
are more typical of the problems encountered in theoret- 
ical physics. 

As an example we consider here a series whose terms 
are products of spherical Bcsscl and Hankcl functions (p. 
435 of Ref. @): 



J 1 (Z)^[7T/(2Z)} 1 / 2 J 1+1/2 (Z), 

^ 1) (z) = [7r/(2z)] 1 /2 jff W /2 ( z) . 



(7.1) 
(7.2) 



Here, the index I is a nonnegative integer, and J2+1/2 

and , 2 are Bessel and Hankel functions, respectively 

(pp. 358 and 360 of Ref. Q]). In physical applications, 
the index I usually finds a natural interpretation as an 
angular momentum quantum number. 

Here, we investigate the following model problem: 



exp( 



-y[l-r}) 



y[l - r] 



1=0 



(2l + l)j l (zry)h' t 1) (iy), (7.3) 
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where y is positive and < r < 1. The spherical Bessel 
and Hankel functions ji{iry) and (iy) can also be ex- 
pressed in terms of modified spherical Bessel functions. 
However, we prefer to retain the given notation because 
there are some inconsistencies in the literature regarding 
the prefactors to be assigned to the modified spherical 
Bessel functions, whereas the spherical Bessel and Han- 
kel functions are defined consistently in most textbooks. 

With the help of known properties of Bessel functions 
and Legendre pol ynom ials it can be shown easily that the 
series expansion (7J3) is a special case of the well known 
addition theorem of the Yukawa potential (p. 107 of Ref. 

H), 



exp(— 7w) 



(rp) 



-1/2 



x J2(2l + l)Pi{coscf>)I l+1/2 ( 1 x)K l+1/2 ( 1 y), (7.4) 
z=o 

where w = (x 2 +y 2 — 2xy cos (j)) 1 ^ 2 , < x < y and 7 > 0. 
This addition theorem is also the Green's function of the 
three-dimensional modified Hclmholtz equation. 

First, we want to analyze the behaviour o f the terms 
of the series on the right-hand side of Eq. ( [7.3| ) if the 
index I becomes large. The leading orders of the asymp- 
totic expansions of the spherical Bessel and Hankel func- 
tions ji(iry) and h^ 1 {iy) as I — > 00 with r and y fixed 
can be obtained easily from the series expansion for ji 
(Eq. (10.1.2) of Ref. || or Eq. (E.17) of Ref. @) and 

from the explicit expression for (Eqs. (10.1.3) and 
(10.1.16) of Ref. |H or Eq. (E.18) of Ref. ||), yielding 

J'(^) ~ ( 2 ?+ 1)!! [1 + Q(1/ ^ )] ' (7 ' 5) 
hf\iy) ~ -i { [l + O (1/0] . (7.6) 

Thus, we obtain for the leading order of the asymptotic 
expansion (I — > 00) of the product of the spherical Bessel 
and the spherical Hankel function: 



ji(iry) hf\iy) 



r l e -y 
(2/ + 1)2/ 



[1 + 0(1/0]. (7.7) 



This asymptotic result sh ows that the infinite series on 
the right-hand side of Eq. ( |7.3| ) converges only for \r\ < 1, 
and that the function exp(— y[l — r])/(y[l — r]) has a 
singularity at r = 1. 

Thus, it is interesting to use the CNCT for an evalu- 
ation of the series on the right-hand side of Eq. ( ff^ ) in 
the immediate vicinity of the singularity. This series is 
similar to the sum over angular momenta encountered in 
QED bound state calculations (cf. Ref. and Eq. (4.2) 
in Ref. @). 

In principle, the spherical Bessel and Hankel functions 
in the products ji(iry) h^\iy) can be evaluated recur- 
sively. If, however, the CNCT is used, this is not useful. 



The products belonging to all angular momenta / are not 
needed. Instead, only the products for some specific an- 
gular momenta have to be evaluated. Accordingly, we 
prefer to evaluate the spherical Bessel function ji(iry) 
from its series expansion (Eq. (10.1.2) of Ref. |34|] or Eq. 
(E.17) of Ref. |8j||) and the spherical Hankel function 

h, {iy) via its explicit expression (Eq. (10.1.16) of Ref. 
_ or Eq. (E.18) of Ref. @). 

It should be taken into account that the recursive eval- 
uation of the terms of angular momentum decomposi- 
tions is normally possible only in the case of simple model 
problems. We shortly outline the applications of the 
CNCT to a more involved computation al p roblem oc- 
curring in QED bound state calculations |?0| . The series 
to be evaluated in the QED calculation is given by Eq. 
(4.2) in Ref. |||. This series can be rewritten as 

00 

S{r,y,t,j) = J(r,y,t,j)^2T K (r,y,t,j) , (7.8) 

K=l 

where J(r, y, t, 7) has a simple mathematical structure. 
The term T K (r, y, t, 7) (k = 1,2,...) can be rewritten as 
(see Eq. (5.7) in Ref. §1) 



xG 



where 



ryy i ( 1 



r y y i ( \ 



a a 2 V t 



\a/ V a a/ V a / 

\a/ V a a/ 



a a 2 V t 



• the functions fi are the radial Dirac-Coulomb wave 
functions (see Eq. (A. 8) in Ref. |PJ), 

• the functions Gg are related to the radial compo- 
nents of the relativistic Green function of the bound 
electron (see Eq. (5.5) in Ref. [§), and 

• the functions A K and A l J are related to the angular 
momentum decomposition of the Green function of 
the virtual photon and are defined in Eq. (5.8) of 
Ref. @. 

The quantity a is a scaling variable for the subsequent 
radial integration over y (see Eq. (4.1) in Ref. Q). The 
variable r denotes the ratio of the two radial coordinates 
(0 < r < 1), and t is related to the (complex) energy 
of the virtual photon (0 < t < 1). The dependence on 
the coupling 7 of the electron to the c entr al field, which 
appears on the left-hand side of Eq. ( |7.9| ), is implicitly 
contained in the functions Gg K on the right-hand side of 
Eq. ([7\|) (see Eq. (A.16) in Ref. @). The propagator 
Gb in Eq. (|7.9|) contains the relativistic Dirac-Coulomb 
Green function, whose radial components are given in 
Eq. (A.16) in Ref. Recurrence formulas relating the 



11 



Gg K for different values of k are not known. 



Therefore 

the terms in the series (7.8) cannot be evaluated by re- 
cursion, but each one of these terms can be computed 
using techniques described in Refs. |30|,|3|]. Note that 
the angular momentum decomposition of Eq. (7.8) docs 



not correspond to an expansion in terms of the QED per- 
turbation theory. The perturbative parameter in QED is 
the elementary charge e or the fine structure constant^] 



a = e 2 /(4 7r). The series (7.8) occurs in the evaluation 



of an energy shift (the Lamb shift) of atomic levels due 
to the self energy of the electron. The self energy is 
an effect described by second-order perturbation theory 
within the framework of QE D. The asymptotic behaviour 
of the terms T K in Eq. (|7.i 




(7.10) 



is similar to that of the model series ( |7.7| ) a nd t o that of 
the Lii-series discussed in Sec. |v| (see Eq. (5.5) and the 
results presented in Table |v|). Slow convergence of the 



series m 



Eq. (|7 



is observed in the limit 



1. 



Returning to our model problem given in Eq. (7.3), we 
consider here the example of r = 0.9999 and y = 0.7. 
This yields 



10~ 5 ^2(21 + 1) 0.9999 x 0.7) h] 1 ' (i 0.7) 



;=o 



-0.142 847 143 207 135. 



(7.11) 



Evaluation of the series over the spherical Bessel and 
Han kel functions via the CNCT is presented in Table 
XII. If this series is evaluated by adding up its terms, 



then about 450 000 products of Bessel and Hankel func- 
tions have to be evaluated for a relative accuracy of 10~ 16 
in the final result. This is to be compared with the ap- 
proximately 300 evaluations of products of Bessel and 
Hankel functions which are needed to compute the 25 
Van Wijngaarden terms presented in Table XII with the 
necessary accuracy. According to our experience, the use 
of the CNCT in QED bound state calculations reduces 
the amount of computer time for slowly convergent series 
of the type of Eq. (|7.8|) by three orders of magnitude. 



Insert Table LXiJJ here 



VIII. SUMMARY AND CONCLUSIONS 

Our approach for the acceleration of the convergence 
of monotone series consists of two steps. In the first step, 



We use natural units as it is customary for QED bound 
state calculations (K = c = eo = 1). 



t he V an Wijngaarden condensation transformation, Eqs. 
( |2.2|) - (0), is applied to the original monotone series. 
In the second step , the convergence of the resulting alter- 
nating series (|T|) is accelerated with the help of suitable 
nonlinear sequence transformations. 

Dan iel [pi[ showed that the transformed alternating 
series ( |2.3| ) converges under relatively mild conditions to 
the same limit as the original series. The condensation 
transformation cannot solve the problems due to slow 
convergence on its own since the transformed alternating 
series does n ot co nverge more rapidly than the original se- 
ries (see Eq. ( |3.2D ). However, the convergence of alternat- 
ing series can be accelerated much more effectively than 
the convergence of the original monotone series. More- 
over, the transformation of alternating series by sequence 
transformations is numerically a relatively stable process, 
whereas the direct transformation of monotone series in- 
evitably leads to the loss of accuracy due to round-off. 

Many sequence transformations are known which are 
in principle capable of accelerating the convergence of al- 
ternating series |j]J|||- However, their efficiency varies 
considerably. In all numerical examples studied in the 
context of this paper, we found that certain sequence 
transformations ||[l2| (which also use explicit estimates 
for the truncation errors as input data) are significantly 
more powerful than other, better known sequence trans- 
formations as for instance Wynn's epsilon algorithm |ll]] 
or the Euler transformation |3{| (which only use the 
partial sums of the infinite series to be transformed 
as input data). Consequently, we use in this article 
only the sequence transformations £l (/?, s„, w„) and 
uj n ) defined in Eqs. (|3~8|) and (EL^]in combi- 



nation with the truncation error estimate 




Those 



properties of these sequence transformations, which are 
needed to apply them successfully in a convergence accel- 
eration process, are described in Section III. In Appendix 
[A|, we explain why these transformations are much more 
effective in the case of alternating series than in the case 
of monotone series, and in Appendix |b| we discuss exact- 
ness properties of these transformations. 

We consistently observed that the combined transfor- 
mation is a remarkably stable numerical process and en- 
tails virtually no loss of numerical significance at interme- 
diate stages. The evaluation of the condensed series (2.4) 
is a computationally simple task, provided that the eval- 
uation of terms of sufficiently high indices in the original 
series is feasible. The nonlinear sequence transformations 
used in the second step are also remarkably stable. The 
final results, which we present in the tables, indicate that 
a relative accuracy of 10~ 14 can be achieved with float- 
ing point arithmetic of 16 decimal digits. Moreover, our 
method remains stable in higher transformation orders. 

In Sections IV - VI we demonstrate that our approach 
is very useful for the computation of special functions 
which are defined by logarithmically convergent series. 
Our approach also works very well if the argument of a 
power series with monotone coefficients is very close to 
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the boundary of the circle of convergence (even in the 
vicinity of singularities). 

In Section IV, we treat the Dirichlet series dlH ) for the 
Riemann zeta function which because of its simplicity is 
excellently suited to illustrate the mechanism of our com- 
bined transformation. The application of the Van Wi- 
jngaarden t rans formation to the logarithmically conver- 
gent series ( L4 ) yields the known alternating series ( |4.3| ) 
which does not c onve rge more rapidly than the original 
monotone series (1.4). However, the alternating series 
can be tr ansfo rme d effe ctively by the sequence transfor- 
mations (3.13) or (3.14) in the case of slow convergence 
or even summed in the case of divergence. Consequently, 
the Riemann zeta function can be evaluated effectively 
and reliably even if its argument z is only slight ly la rger 
than one (Table ||) , if both the monotone series ( |l.4[ ) and 
the alternating series (4.3) diverge (Table III), or if the 
argument z is complex (Tabic IV). 

In Section [v| we treat the Lerch transcendent $(z, s, a) 
which contains many other special functions as special 
cases, for example the Riemann zeta function or the 
poly logarithms. In Tables |v| - |VIII we show that the 
combined transformation makes it possible to evaluate 



&(z, s, a) effectively and reliable via the power series (5.1) 



even if z is very close to the boundary of the circle of con- 
vergence. 

In Section [v] we discuss the evaluation of a nontcr- 
minating generalized h ype rgeometric series p +\F p via its 
defining power series (1.3) which converges for \z\ < 1. 
In the case of a Gaussian hypergeometric series %Fi, ex- 
plicit analytic continuation formulas are known which 
transform a hypergeometric series with argument z into 
the sum of two hypergeometric series 2-F1 with argument 
1 — z. Thus, if the convergence of a hypergeometric series 
2-F1 is slow because its argument z is only slightly smaller 
then one, then the two transformed hypergeometric series 
2-F1 with argument 1 — z will converge rapidly. Moreover, 
if a Gaussian hypergeometric series %Fi with unit argu- 
ment z = 1 converges, the n its value is given by the Gauss 
summation theorem (S.5). In the case of a generalized 
hypergeometric series P +iF p with p > 2, the situation is 
much more difficult. Explicit analytic continuation for- 
mulas are either unknown or they become increasingly 
complicated with increasing p. Similarl y, s imple ana- 
logues for the Gauss summation theorem (3.5) exist only 
for certain values of p and for certain co mbi nations of the 
parameters of the p +iF p . In Tables IX - XI we show that 
a generalized hypergeometric series P +iF p (j> > 2) with 
an argument z that is only slightly smaller than one or 
with unit argument [z = 1) can be computed effectively 
and reliably with the help of the combined transforma- 
tion. 

Partial wave decompositions of Green's functions, 
which occur for example in quantum electrodynamic 
bound state calculations, entail more complex mathemat- 
ical entities than series expansions for special functions. 
Nevertheless, we show in Section VII that the combined 



tial wave decompositions. The series over products of 
spherical Bessel and Hankel functions considered in this 
paper serves as a model problem for the angular momen- 
tum decomposition of more complex Green's functions, 
as for example the relativistic Green's function of the 
bound electron [Q. In the model problem studied here 
as well as in a recent evaluation of self energy correc- 
tions in bound systems |5(j, we observed a reduction in 
computer time by three orders of magnitude. Thus, the 
combined transformations makes extensive and highly ac- 
curate calculations feasible in situations that could oth- 
erwise be too time-consuming. 
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APPENDIX A: ON THE EFFICIENCY OF THE 
TRANSFORMATION OF ALTERNATING AND 
MONOTONE SERIES 

In the case of the Levin transformation £j^(/3, s n ,u) n ), 
Eq. ( |3.8| ) , and the closely related Weniger transformation 
uj n ), Eq. (|3~9|), it is relatively easy to under- 
stand why alternating series can be transformed much 
more effectively than monotone series. These two trans- 
formations are both special cases of the following trans- 
formation which is characterized by the positive weights 
w k (n): 



T^ n) (w k (n); s n ,uj n ) 



A k {w k (n)s n /uj n } 
A k {w k (n)/uj n } 



(Al) 



If Wk(n) = (n + 0) we obtain Levin's transforma- 
tion, and it Wk(n) = (n + 0)k-i, we obtain Weniger's 
transformation. 

Since the difference operator A k is linear, Eq. ( |Al[ ) can 
also be rewritten as follows: 



T^ n \w k (n);s ni uJn) = s + 



A k {w k (n)[s n - s]/uj n } 
A k {w k (n)/uj n } 



(A2) 



transformation can be applied successfully to these par- 
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Obviously, the sequence transformation T k converges 
to the (generalized) limit s of the input sequence {s n }^L 
if the remainder estimates {w n }%L can be chosen in such 
a way that the ratio on the right-hand side becomes neg- 
ligibly small. 

Loosely speaking, this means that we have to choose 
the remainder estimates in such a way tha t th e numerator 
of the ratio on the right-hand side of Eq. ( A2) becomes as 
small as possible and the denominator becomes as large 
as possible. 

If we can find remainder estimates such that 



LO n [c+ 0(n *)] , n 



(A3) 



then the weighted difference operator A k Wk(n) will make 
the numerator small, no matter whether the input data 
{s„}^Lq are the partial sums of an alternating or of a 
monotone series. 

The situation is quite different in the case of the de- 
nominator. Application of the weighted difference oper- 
ator A k Wk(n) to 1/wn yields 

A k {w k (n)/uj n } 



(A4) 



If the remainder estimates all have the same sign, the 
denominator will also become relatively small because of 
cancellation due to differencing (hopefully not as small 
as the numerator because then the transformation pro- 
cess would not converge). This cancellation process is 
also the major source of numerical instabilities which are 
magnified since they occur in the denominator. 

If the remainder estimates are strictly alternating in 
sign, io n = (— l) n |w n |, then we obtain (w k (n) is by as- 
sumption positive) 



A k {w k (n)/w n } = (-l) k+n 



k\ w k (n + j) 



Jn+j | 



(A5) 



Thus, all terms of the denominator sum have the same 
sign, and there is no cancellation as in the case of mono- 
tone remainder estimates and also no source of numeri- 
cal instabilities. Consequently, the denominator becomes 
relatively large which improves convergence. 



As an example, we apply the transformation (Al) to 
the monotone model sequence 



Cl 



(n + l) a {n + l) a+1 
as well as to the alternating model sequence 



(A6) 



+ (-!)" 



Co 



Cl 



(n+l) a (n+l) a+1 



(A7) 



and derive asymptotic (n — > 00) transformation error 
estimates. 



For that purpose, we assume a > - which implies 
that the two model sequences converge to their limits s 
and t, respectively - and w k (n) = 0(n fc_1 ) as n — > 00. 
In the case of the monotone input sequence {sn}^Lo> we 
choose u> n = (n+ 1)~ Q , and in the case of the alternating 
input sequence {t n }^ =Q , we choose ui n = (— l) n (n+l)~ a . 
Then, we obtain for the numerators 



k ,w k (n)[s n - s] 



= A 



(n + 1)" Q 

k ( w k (n)[t n -t] 



(-l)™(n + 1)~ Q 



j = 0(n" fe - 1 ). (A8) 



The estimates for the denominators differ cons ider ably. 
In the case of the monotone input sequence (A6), we 
obtain 



A 



w k {n) 



(n+iy 



0(n a ~ 



(A9) 



which in combination with Eq. ( |A8| ) yields the following 
asymptotic transformation error estimate as n — > 00: 



T k {n) (w k (n);s n ,(n+l) a ) - s = Q ( n -k\ 



(A10) 



In the case of the alternating input sequence (A7), we 
exploit the fact that the first term with j = in the sum 
on the right-hand side of Eq. ( A5) is smaller in magnitude 
than the whole sum, yielding 



\A k {w k (n) /[(-!}" (n + l)-<*}}\ 
1 



< 



\w k (n)/[(-l) n (n+ l)~ a }\ ' 



(All) 



Of course, we could also try to construct more sophis- 
ticated estimates for the denominators. However, the 



relatively crude estimate (All) suffices f or o ur purposes 
since it implies in combination with Eq. (A8) 



T k {n \w k (n);t n ,(-iy(n 



t n t 



0{n- 2k ) 



(A12) 



The two transformation error estimates ( A10| ) and ( A12| ) , 
which both hold as n — * 00, show that there is a substan- 
tial difference between the transformation of monotone 
and alternating series. There is a considerable amount of 
evidence that this conclusion is actually generally true, 
i.e., also in the case of Pade approximants and other se- 
quence transformations. 



APPENDIX B: EXACTNESS RESULTS 



In this Appendix, we analyze exactness properties 
of the Levin transformation [j3, s n , uj n ), Eq. (|3.g|), 
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and of the Weniger transformation Sj^{j3,Sn,u} n ), Eq. 



(3.9), and their variants d k n '((3, S n ), Eq. ( 3.13| ), and 



Sr'(p,S n ), Eq. j04| ), which both use the remainder 



estimate ( 3.12 ) proposed by Smith and Ford 

For that purpose we introduce in Eq. (A2) the remain- 
der r n = s n — s and obtain 



(w k {n); s n ,oj n 



A k {wk{n)r n /uj n } 
A k {w k (n)/uj n } 



(Bl) 



Obviously, the general sequence transformation 77 , Eq. 



(Al), which contains C k n \f3, s n ,u) n ) and 5^ (/?, s n ,u> n ) 
as special cases, is exact for a given input se- 
quence {s„}^L ^ the difference operator A fc annihilates 
w k (n)r n /u> n but not w k (n)/uj n . 



Thus, if we want to prove the exactness of T k 



(n) 



for 



some sequence {s„}^L , we need to know explicit expres- 
sions for the remainders {rn}n°=o °f the m put sequence 
as functions of n. 

In the case of the alternating series (flj]) for the Rie- 
mann zeta function with zero or negative integral argu- 
ment, t his can be accomplished easily. Replacing z by — I 
in Eq. (O) with I = 0, 1, 2, ... we obtain 



CO 



1 oo 

-0 = t^ti£(-W + 1 )' 



(B2) 



For 1 = 0, this series is the negative of the geometric 
series 1/(1 + x) = Y^^f{~ x Y at a; = 1. According to 
Theorem 12-9 of Ref. M, the sequence transformations 
C k n \f3, s n ,uj n ) and sj^{/3,s n ,u n ) are exact for the geo- 
metric series if the first term neglected in the partial sum 
is used as remainder estim ate (which corresponds to the 
remainder estimate ( 3.12 )). Thus, the exactness of these 
sequence transformations for the infinite series ( |B2| ) has 
to be analyzed only for I > 1. 

For the determination of its truncation error with ar- 
bitrary integral I > 1, we rewrite the infinite series ( ]b"2| ) 
as follows: 



3=0 



= £(-iyc7+iy + £ c-iyci+i) 

n 

= ^(-iy(j+iy 



(B3) 



3=0 



{ _ ir+ isr { _ ir{n+u+2) i 



(B4) 



The infinite series on the right-hand side of Eq. (B4) 
obviously diverges. However, it can be summed easily 
since it can be represented as a derivative of the geometric 
series 1/(1 + x) = J2jLo(^ x y at x = 1: 



]T(-l)> + ^ + 2)< 

{oo 
V (-l) u {n + v + 2) l x 
Z — ✓ 




(B5) 
(B6) 
(B7) 



By inserting this into Eq. (B2) we obtain the following 
representation of £(— I) as a partial sum plus an explicit 
expression for the truncation error: 



1 f - 

CH) = yt^tt EH^ + i)' 



+ (-1) 



n+l 



d V X n+1 



dx 



l + x 



(B8) 



If we set in Eq. (B8) n = —1, then the partial sum is an 
empty sum and vanishes, yielding the following explicit 
expression for the Ricmann zeta function with zero or 
negative integer argument: 



CH) 



i 



1 - 2'+! 



d 

dx 



1 



l+x 



(B9) 



Equivalent explicit expressions can be obtained via the 
substitution y = ln(x). Obviously, we have 



dy + 1 



y=\a(x) 



(BIO) 



In this way, we obtain: 

C(-0 



1 



1 - 2'+! 



1 1 



dy 



1 



1 + e^ 



(Bll) 



1 f - 



y=o 



+ (-1) 



n+l 



d y e (" +i ^ 

dy + / 1 + e^ 



(B12) 



Combination of Eqs. ( |4.11b| ), (|B9|), and ( |Bll| ) yields the 
following expressions for the Bernoulli numbers with even 
indices: 



B21 = 



-21 
1 - 2 21 

-21 



d V 1 - 1 1 



dx 
d 



l + x 

21-1 



1 



i + cy 



(B13) 
(B14) 



y=o 
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It follows from Eq. (Bl) that the general sequence 



transformation T k 



(n) 



Eq. (Al), is exact for some in- 
put sequence {s n }n°=o if the expression Wk(n)r n /u n with 
T n = s n — s is a polynomial of degree k — 1 in n. Thus, 
Eq. flB7|) implies that in the case of the Weniger transfor- 
mation S n (l, Sp) , Eq. ( 3.14 ), which uses the remainder 
estimate (3.12), we have to show that the ratio 



i) fc -iIT=o(-i)> 



2 y 



(-1)"+! {n + 2) 1 



(B15) 



is for sufficiently large values of k a polynomial of degree 
k - 1 in ra. For I = 1, 2, 3, and 4, Eq. (|B7|) yields the 
following explicit expressions: 



^(-inn+, + 2) 1 

oo 

]T(-l)> + ^ + 2) 2 



2n + 3 



(n +!)(» + 2) 



(B16) 



, (B17) 



u=0 



£(-!)> + v + 2f 



u=0 



(2n + 3)(2n 2 + 6n + 3) 



£(-!)> + „ + 2) 



(n + l)(n + 2)(n 2 + 3n+l) 



(B18) 



(B19) 



Since (n + = (n + l)(n + 2) . . . (n + k — 1), the ratio 
(B15) is for I = 1, 2 and k > 3 a polynomial of degree 
fe — 1 in n, whereas for I — 3, 4 it is a rational expression 
in n which is not annihilated by the operator A k . 

Thus, the exactness of the transformation 5^(1, So) 
for £(— 1) as in Table III, or for £(— 2) is more or less 



accidental. If, however, we choose in the Levin transfor- 
mation dn^ (J3, So), Eq. ( 3.13 ), with (3 = 2 instead of the 
usual (3 — 1, we find that the corresponding ratio 



(n+2) fc - i Er = o(-ir( 



2)' 



(_l)n+l ( n + 2) ; 



(B20) 



is a polynomial of degree k — 1 in n for all Z = 0, 1, 2, ... 
and fc > l + l. This follows at once from the fact that the 



differential operator in Eq. (B7) produces a polynomial 
of degree ! inn. 

In the same way, it can be shown that the Weniger 
transformation S„ (/3, So), Eq. ( |3.14 ), with (3 = 2 is ex- 
act for the hypergeometric series for all I = 0, 1, 2, . . . 
and for k > I + 1: 



niF (i + l;-l) - ^(-l)^(j + l) 

3=0 

This can also be rewritten as follows: 



2 l + 



i ■ (B21) 



ni.Fo(f + l;-l) 
= £(-1)^ + 1)* 

OO 

+ ( _ ir+ i^ ( _ ir(r 



■2)j 



(B22) 



(-1)^(7 + 



+ (-1)" +1 (n + 2)j 2 Fi(l, n + i + 2; n + 2; -1) (B23) 



3=0 



+ (-l) n+1 E H)i(n + J +2)w(l/2) 

In this case, the corresponding ratio 

("+2)t-iEM(-lf(" + ^2)i 

(_l)n+l ( n + 2 ) z 



(B24) 



(B25) 



is a polynomial of degree fc — 1 in n for alH = 0, 1, 2, ... 
and for k > l + l. This follows at once from the fact that 



the second sum on the right-hand side of Eq. (B24) is a 
polynomial of degree I in n. 



However, the value f3 = 2 in the Eqs. (3.13) and (3.14) 
cannot significantly improve convergence for the zeta 
function. We also investigated the performance of the 
transformations dn°' ) (2,S„) and <5i 0) (2,S rl ). Except for 
the cases of z = — 1 or z = — 2 in which the Levin trans- 
formation becomes exact with (3 = 2, the transformations 
dn^ (2, S n ) and Sn^ (2, S„) yield results which differ only 
marginally from the results obtained by dn^ (l, S n ) and 
5^ (l, S„), which are presented in Tables || - |fv| . 
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TABLES 



TABLE I. Demonstration that Van Wijngaarden's transformation, Eqs. ( p.2[ ) - (2.4), is a rear- 
rangement of the original series y"]^_ n a(k). 



(-l)'A, 






Leading terms of the original 


series 




A 


o(0) 


+2a(l) 




+4a(3) 




+ 


-Ar 




-a(l) 




-2 o(3) 






A 2 






+o(2) 






+ 


-A 3 








-o(3) 






A 4 










+a(4) 






a(0) 


+a(l) 


+o(2) 


+a(3) 


+o(4) 


+ 



TA BLE II. Evaluation of the Riemann zeta function of argument 2 = 1.01 by accelerating the convergence of the alternating 
series (4.3). The result is given as io _3 c(i.oi). 



n 








3» 










E n 










dl 0) ( 


i,s 


0) 








5i 0) ( 


1,S ) 










144 


770 


081 


711 


084 





144 


770 


081 


711 


084 





144 


770 


081 


711 


084 





144 


770 


081 


711 


084 


1 





072 


885 


040 


855 


542 





090 


606 


301 


069 


428 





101 


569 


133 


143 


252 





101 


569 


133 


143 


252 


2 





120 


614 


482 


322 


980 





096 


697 


481 


252 


857 





100 


456 


642 


533 


059 





100 


456 


642 


533 


059 


3 





084 


920 


235 


019 


068 





098 


985 


546 


018 


036 





100 


587 


783 


459 


042 





100 


579 


332 


613 


649 


4 





113 


411 


984 


373 


518 





099 


901 


837 


494 


047 





100 


577 


428 


866 


203 





100 


578 


083 


572 


921 


5 





089 


712 


109 


307 


161 





100 


283 


957 


662 


399 





100 


577 


954 


415 


585 





100 


577 


949 


566 


834 


6 





109 


994 


997 


614 


328 





100 


447 


835 


715 


031 





100 


577 


944 


116 


204 





100 


577 


943 


567 


122 


7 





092 


271 


153 


050 


434 





100 


519 


572 


572 


454 





100 


577 


943 


249 


050 





100 


577 


943 


346 


150 


8 





108 


007 


136 


313 


467 





100 


551 


470 


653 


282 





100 


577 


943 


342 


049 





100 


577 


943 


338 


734 


9 





093 


859 


665 


080 


579 





100 


565 


830 


599 


811 





100 


577 


943 


338 


553 





100 


577 


943 


338 


503 


10 





106 


708 


750 


240 


925 





100 


572 


360 


153 


981 





100 


577 


943 


338 


482 





100 


577 


943 


338 


497 


11 





094 


940 


666 


205 


327 





100 


575 


353 


801 


276 





100 


577 


943 


338 


498 





100 


577 


943 


338 


497 


12 





105 


794 


821 


569 


601 





100 


576 


735 


870 


616 





100 


577 


943 


338 


497 





100 


577 


943 


338 


497 


13 





095 


723 


429 


487 


784 





100 


577 


377 


707 


430 





100 


577 


943 


338 


497 





100 


577 


943 


338 


497 


14 





105 


116 


912 


361 


076 





100 


577 


677 


299 


954 





100 


577 


943 


338 


497 





100 


577 


943 


338 


497 


15 





096 


316 


203 


847 


728 





100 


577 


817 


763 


434 





100 


577 


943 


338 


497 





100 


577 


943 


338 


497 



exact 



0.100 577 943 338 497 



0.100 577 943 338 497 



0.100 577 943 338 497 



0.100 577 943 338 497 
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TABLE III. Evaluation of the Riemann zeta function of argument z = —1 by summing the 
divergent alternating series (4. J). The result is given as 10£(— 1). 



n 






S n 








a 


ffll / 
n \ 


1,S 


) 






5 


i 0) (l,S ) 







-3 


333 


333 


333 


333 


333 


-3 


333 


333 


333 


333 


333 


-3 


333 


333 


333 


333 


333 


1 


3 


333 


333 


333 


333 


333 


-0 


666 


666 


666 


666 


667 


-0 


666 


666 


666 


666 


667 


2 


-6 


666 


666 


666 


666 


667 


-0 


860 


215 


053 


763 


441 


-0 


860 


215 


053 


763 


441 


3 


6 


666 


666 


666 


666 


667 


-0 


830 


449 


826 


989 


619 


-0 


833 


333 


333 


333 


333 


4 


-10 


000 


000 


000 


000 


000 


-0 


833 


557 


890 


954 


819 


-0 


833 


333 


333 


333 


333 


5 


10 


000 


000 


000 


000 


000 


-0 


833 


319 


627 


418 


409 


-0 


833 


333 


333 


333 


333 


6 


-13 


333 


333 


333 


333 


333 


-0 


833 


334 


020 


666 


741 


-0 


833 


333 


333 


333 


333 


7 


13 


333 


333 


333 


333 


333 


-0 


833 


333 


304 


093 


535 


-0 


833 


333 


333 


333 


333 


8 


-16 


666 


666 


666 


666 


667 


-0 


833 


333 


334 


413 


139 


-0 


833 


333 


333 


333 


333 


9 


16 


666 


666 


666 


666 


667 


— 


833 


333 


333 


298 


109 


— 


833 


333 


333 


333 


333 


10 


-20 


000 


000 


000 


000 


000 


-0 


833 


333 


333 


334 


362 


-0 


833 


333 


333 


333 


333 


11 


20 


000 


000 


000 


000 


000 


-0 


833 


333 


333 


333 


306 


-0 


833 


333 


333 


333 


333 


12 


-23 


333 


333 


333 


333 


333 


-0 


833 


333 


333 


333 


334 


-0 


833 


333 


333 


333 


333 


13 


23 


333 


333 


333 


333 


330 


-0 


833 


333 


333 


333 


333 


-0 


833 


333 


333 


333 


333 


14 


-26 


666 


666 


666 


666 


670 


-0 


833 


333 


333 


333 


333 


-0 


833 


333 


333 


333 


333 


15 


26 


666 


666 


666 


666 


670 


-0 


833 


333 


333 


333 


333 


-0 


833 


333 


333 


333 


333 



exact -0.833 333 333 333 333 -0.833 333 333 333 333 -0.833 333 333 333 333 



TABLE IV. Evaluation of C(l/2 + 13 -7 i) with the CNCT. 



n 












d {0) 


(l,S ) 
























<5i 0) 






















414 


107 


543 


949 


134 


+ 





017 


316 


297 


125 


790 


i 





414 


107 


543 


949 


134 


+ 





017 


316 


297 


125 


790 i 


1 





575 


871 


239 


097 


112 







042 


690 


435 


565 


758 


i 





575 


871 


239 


097 


112 







042 


690 


435 


565 


758 i 


2 





523 


424 


912 


174 


020 


+ 





152 


835 


043 


959 


961 


i 





567 


958 


887 


269 


553 


+ 





129 


913 


386 


486 


220 i 


3 





481 


953 


715 


196 


159 







288 


400 


086 


031 


923 


i 





474 


129 


917 


411 


618 







168 


775 


423 


138 


291 i 


4 





012 


442 


899 


246 


184 







237 


603 


260 


694 


125 


i 


-0 


180 


827 


868 


994 


142 







367 


542 


940 


737 


051 i 


5 





123 


074 


021 


609 


358 







316 


357 


718 


264 


423 


i 





126 


392 


529 


409 


594 







290 


235 


127 


404 


228 i 


6 





105 


377 


569 


236 


175 







313 


246 


538 


788 


829 


i 





107 


386 


234 


787 


298 







317 


087 


400 


005 


856 i 


7 





107 


635 


288 


132 


001 







312 


843 


180 


304 


712 


i 





107 


124 


241 


668 


490 







312 


622 


662 


817 


549 i 


8 





107 


429 


326 


957 


578 







312 


999 


708 


812 


577 







107 


489 


121 


873 


498 







312 


978 


336 


617 


038 i 


9 





107 


438 


933 


679 


469 







312 


974 


184 


308 


675 


i 





107 


436 


183 


812 


222 







312 


980 


180 


873 


835 i 


10 





107 


439 


640 


888 


613 







312 


976 


813 


188 


762 


i 





107 


439 


393 


557 


558 







312 


976 


229 


866 


877 i 


11 





107 


439 


434 


156 


190 







312 


976 


661 


398 


796 


i 





107 


439 


488 


450 


447 







312 


976 


678 


675 


422 i 


12 





107 


439 


457 


152 


512 







312 


976 


659 


181 


027 


i 





107 


439 


453 


065 


793 







312 


976 


661 


841 


904 i 


13 





107 


439 


455 


840 


151 







312 


976 


660 


717 


255 


i 





107 


439 


455 


877 


677 







312 


976 


660 


319 


609 i 


14 





107 


439 


455 


825 


365 







312 


976 


660 


547 


552 


i 





107 


439 


455 


848 


368 







312 


976 


660 


568 


977 % 


15 





107 


439 


455 


836 


355 







312 


976 


660 


556 


014 


i 





107 


439 


455 


833 


989 







312 


976 


660 


556 


440 i 


16 





107 


439 


455 


835 


269 







312 


976 


660 


556 


233 


i 





107 


439 


455 


835 


348 







312 


976 


660 


556 


072 i 


17 





107 


439 


455 


835 


311 







312 


976 


660 


556 


157 


i 





107 


439 


455 


835 


317 







312 


976 


660 


556 


169 i 


18 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 


i 





107 


439 


455 


835 


312 







312 


976 


660 


556 


163 i 


19 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 


i 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 i 


20 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 


i 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 i 


21 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 


i 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 i 


22 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 


i 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 i 


23 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 


i 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 i 


24 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 


i 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 i 


25 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 


i 





107 


439 


455 


835 


313 







312 


976 


660 


556 


163 i 


exact 


0.107 439 455 835 313 




0.312 976 660 556 163 


i 


0.107 439 455 835 313 




0.312 976 660 556 163 i 
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TABLE V. Evaluation of 1CT 2 Lii (0.99999) = -1(T 2 ln(O.OOOOl) with the CNCT. 



n 






1 


3„ 










di 0) ( 


>,S 


o) 








5i 0) ( 


l,So) 







0. 


162 


768 


973 


713 


089 


0. 


.162 


768 


973 


713 


089 


0. 


162 


768 


973 


713 


089 


1 


0. 


086 


384 


436 


856 


544 


0. 


.116 


225 


388 


785 


336 


0. 


116 


225 


388 


785 


336 


2 


0. 


135 


357 


615 


351 


010 


0. 


.114 


995 


657 


006 


664 


0. 


114 


995 


657 


006 


664 


3 


0. 


099 


665 


296 


922 


988 


0. 


.115 


140 


148 


148 


939 


0. 


115 


131 


002 


772 


470 


4 


0. 


127 


575 


317 


431 


702 


0. 


.115 


128 


665 


188 


679 


0. 


115 


129 


400 


970 


919 


5 


0. 


104 


755 


344 


851 


635 


0. 


.115 


129 


271 


216 


942 


0. 


115 


129 


261 


517 


373 


6 


0. 


123 


997 


631 


730 


577 


0. 


.115 


129 


254 


851 


082 


0. 


115 


129 


254 


830 


226 


7 


0. 


107 


401 


422 


517 


315 


0. 


.115 


129 


254 


746 


908 


0. 


115 


129 


254 


725 


403 


8 


0. 


121 


964 


815 


378 


398 


0. 


.115 


129 


254 


612 


355 


0. 


115 


129 


254 


664 


668 


9 


0. 


109 


009 


755 


125 


042 


0. 


.115 


129 


254 


647 


759 


0. 


115 


129 


254 


647 


723 


10 


0. 


120 


662 


087 


466 


136 


0. 


.115 


129 


254 


650 


661 


0. 


115 


129 


254 


648 


068 


11 


0. 


110 


085 


384 


510 


686 


0. 


.115 


129 


254 


650 


507 


0. 


115 


129 


254 


649 


602 


12 


0. 


119 


759 


673 


589 


890 


0. 


.115 


129 


254 


649 


358 


0. 


115 


129 


254 


649 


772 


13 


0. 


110 


852 


765 


866 


205 


0. 


.115 


129 


254 


649 


714 


0. 


115 


129 


OK A 




711 


14 


0. 


119 


099 


529 


780 


039 


0. 


.115 


129 


254 


649 


727 


0. 


115 


129 


254 


649 


700 


15 


0. 


111 


426 


375 


175 


159 


0. 


.115 


129 


254 


649 


696 


0. 


115 


129 


254 


649 


702 


16 


0. 


118 


596 


725 


346 


641 


0. 


.115 


129 


254 


649 


702 


0. 


115 


129 


254 


649 


702 


17 


0. 


111 


870 


534 


473 


655 


0. 


.115 


129 


254 


649 


703 


0. 


115 


129 


254 


649 


702 


18 


0. 


118 


201 


666 


639 


655 


0. 


.115 


129 


254 


649 


702 


0. 


115 


129 


254 


649 


702 


19 


0. 


112 


224 


086 


515 


227 


0. 


.115 


129 


254 


649 


702 


0. 


115 


129 


254 


649 


702 


20 


0. 


117 


883 


505 


752 


715 


0. 


.115 


129 


254 


649 


702 


0. 


115 


129 


254 


649 


702 



exact 0.115 129 254 649 702 0.115 129 254 649 702 0.115 129 254 649 702 



TABLE VI. Evaluation of 10" 1 Li 2 (0.99999) with the CNCT. 



n 






t 


3„ 












'l,S 


o) 








5i 0) ( 


l,So) 







0. 


199 


982 


280 


324 


442 


0. 


.199 


982 


280 


324 


442 


0. 


199 


982 


280 


324 


442 


1 


0. 


149 


990 


640 


162 


221 


0. 


.165 


371 


886 


328 


955 


0. 


165 


371 


886 


328 


955 


2 


0. 


172 


207 


484 


145 


972 


0. 


.164 


381 


453 


915 


497 


0. 


164 


381 


453 


915 


497 


3 


0. 


159 


711 


414 


066 


111 


0. 


.164 


488 


426 


505 


073 


0. 


164 


482 


760 


527 


739 


4 


0. 


167 


708 


334 


514 


884 


0. 


.164 


480 


538 


599 


000 


0. 


164 


481 


025 


806 


042 


5 


0. 


162 


155 


301 


413 


564 


0. 


.164 


480 


897 


128 


353 


0. 


164 


480 


896 


552 


227 


6 


0. 


166 


234 


803 


732 


817 


0. 


.164 


480 


894 


717 


325 


0. 


164 


480 


893 


640 


937 


7 


0. 


163 


111 


643 


694 


762 


0. 


.164 


480 


893 


606 


728 


0. 


164 


480 


893 


688 


676 


8 


0. 


165 


579 


162 


855 


583 


0. 


.164 


480 


893 


702 


656 


0. 


164 


480 


893 


698 


442 


9 


0. 


163 


580 


602 


633 


196 


0. 


.164 


480 


893 


699 


380 


0. 


164 


480 


893 


699 


234 


10 


0. 


165 


232 


198 


806 


363 


0. 


.164 


480 


893 


699 


272 


0. 


164 


480 


893 


699 


288 


11 


0. 


163 


844 


487 


813 


342 


0. 


.164 


480 


893 


699 


295 


0. 


164 


480 


893 


699 


292 


12 


0. 


165 


026 


841 


360 


043 


0. 


.164 


480 


893 


699 


292 


0. 


164 


480 


893 


699 


293 


13 


0. 


164 


007 


426 


937 


253 


0. 


.164 


480 


893 


699 


293 


0. 


164 


480 


893 


699 


293 


14 


0. 


164 


895 


394 


970 


447 


0. 


.164 


480 


893 


699 


293 


0. 


164 


480 


893 


699 


293 


15 


0. 


164 


115 


002 


453 


607 


0. 


.164 


480 


893 


699 


293 


0. 


164 


480 


893 


699 


293 



exact 0.164 480 893 699 293 0.164 480 893 699 293 0.164 480 893 699 293 
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TABLE VII. Evaluation of 1CT 1 Li 3 (0.99999) with the CNCT. 



n 






1 


3„ 










di 0) ( 


1,S 


") 








5i 0) ( 


1,S ) 







0. 


133 


331 


333 


415 


539 


0. 


.133 


331 


333 


415 


539 


0. 


133 


331 


333 


415 


539 


1 


0. 


116 


665 


166 


707 


769 


0. 


.120 


474 


532 


168 


000 


0. 


120 


474 


532 


168 


000 


2 


0. 


121 


603 


216 


117 


468 


0. 


.120 


176 


326 


936 


846 


0. 


120 


176 


326 


936 


846 


3 


0. 


119 


520 


007 


764 


208 


0. 


.120 


206 


042 


152 


677 


0. 


120 


204 


748 


497 


388 


4 


0. 


120 


586 


594 


446 


594 


0. 


.120 


203 


955 


328 


380 


0. 


120 


204 


079 


128 


106 


5 


0. 


119 


969 


366 


038 


597 


0. 


.120 


204 


046 


033 


711 


0. 


120 


204 


045 


387 


208 


6 


0. 


120 


358 


052 


152 


572 


0. 


.120 


204 


045 


725 


809 


0. 


120 


204 


045 


378 


284 


7 


0. 


120 


097 


666 


726 


411 


0. 


.120 


204 


045 


413 


120 


0. 


120 


204 


045 


434 


802 


8 


0. 


120 


280 


540 


991 


745 


0. 


.120 


204 


045 


439 


707 


0. 


120 


204 


U40 




553 


9 


0. 


120 


147 


227 


650 


952 


0. 


.120 


204 


045 


438 


749 


0. 


120 


204 


045 


438 


726 


10 


0. 


120 


247 


386 


435 


540 


0. 


.120 


204 


045 


438 


729 


0. 


120 


204 


045 


438 


733 


11 


0. 


120 


170 


239 


824 


481 


0. 


.120 


204 


045 


438 


733 


0. 


120 


204 


045 


438 


733 


12 


0. 


120 


230 


916 


813 


857 


0. 


.120 


204 


045 


438 


733 


0. 


120 


204 


045 


438 


733 


13 


0. 


120 


182 


336 


147 


846 


0. 


.120 


204 


045 


438 


733 


0. 


120 


204 


045 


438 


733 


14 


0. 


120 


221 


833 


436 


597 


0. 


.120 


204 


045 


438 


733 


0. 


120 


204 


045 


438 


733 


15 


0. 


120 


189 


289 


161 


290 


0. 


.120 


204 


045 


438 


733 


0. 


120 


204 


045 


438 


733 



exact 0.120 204 045 438 733 0.120 204 045 438 733 0.120 204 045 438 733 



TABLE VIII. Evaluation of 10 4 $(0.99999, 2, 10000) with the CNCT. 



n 






1 


3„ 










di 0) ( 


'l,S 


o) 








5i 0) ( 


l,So) 







1 


.152 


086 


970 


131 


424 


1. 


.152 


086 


970 


131 


424 


1 


,152 


086 


970 


131 


424 


1 





.576 


093 


485 


065 


712 


0. 


806 


478 


876 


912 


452 





,806 


478 


876 


912 


452 


2 


0. 


.960 


055 


803 


546 


361 


0. 


.797 


618 


192 


129 


198 





,797 


618 


192 


129 


198 


3 


0. 


.672 


109 


050 


515 


104 


0. 


.798 


663 


645 


011 


412 





,798 


596 


144 


946 


064 


4 





.902 


446 


455 


721 


367 


0. 


.798 


581 


028 


864 


897 





,798 


586 


253 


716 


867 


5 





.710 


515 


275 


487 


446 


0. 


.798 


585 


227 


987 


408 





,798 


585 


188 


634 


170 


6 





.875 


013 


443 


138 


958 


0. 


.798 


585 


145 


208 


936 





,798 


585 


140 


857 


888 


7 





.731 


090 


035 


137 


739 


0. 


.798 


585 


138 


553 


527 





,798 


585 


139 


249 


075 


8 





.859 


010 


851 


725 


411 


0. 


.798 


585 


139 


276 


063 





,798 


585 


139 


237 


667 


9 





.743 


892 


107 


147 


896 


0. 


798 


585 


139 


218 


618 





,798 


585 


139 


229 


500 


10 


0. 


.848 


536 


431 


310 


159 


0. 


.798 


585 


139 


219 


617 





,798 


585 


139 


222 


908 


11 





.752 


620 


788 


733 


222 


0. 


.798 


585 


139 


223 


310 





,798 


585 


139 


222 


175 


12 





.841 


150 


625 


351 


432 


0. 


.798 


585 


139 


222 


720 





,798 


585 


139 


222 


491 


13 





.758 


951 


478 


583 


304 


0. 


.798 


585 


139 


222 


444 





,798 


585 


139 


222 


559 


14 





.835 


664 


028 


102 


224 


0. 


.798 


585 


139 


222 


555 





,798 


585 


139 


222 


550 


15 





.763 


752 


250 


680 


046 


0. 


.798 


585 


139 


222 


555 





,798 


585 


139 


222 


548 


16 





.831 


428 


053 


244 


539 


0. 


.798 


585 


139 


222 


546 





,798 


585 


139 


222 


548 


17 





.767 


517 


561 


053 


133 


0. 


.798 


585 


139 


222 


548 





,798 


585 


139 


222 


548 


18 


0. 


.828 


059 


092 


035 


324 


0. 


.798 


585 


139 


222 


548 





,798 


585 


139 


222 


548 


19 





.770 


549 


625 


376 


190 


0. 


.798 


585 


139 


222 


548 





,798 


585 


139 


222 


548 


20 


0. 


.825 


315 


796 


529 


872 


0. 


.798 


585 


139 


222 


548 





,798 


585 


139 


222 


548 



exact 0.798 585 139 222 548 0.798 585 139 222 548 0.798 585 139 222 548 
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TABLE IX. Evaluation of 1CT 4 3 F 2 (l, 3/2, 5; 9/8, 47/8; 0.99999) with the CNCT. 



n 






1 


3„ 










di 0) ( 


>,S 


o) 








5 ( ° ] { 


1,S ) 







0. 


.343 


961 


195 


195 


881 


0. 


.343 


961 


195 


195 


881 


0. 


.343 


961 


195 


195 


881 


1 





.172 


030 


597 


597 


940 


0. 


.240 


789 


533 


227 


428 





,240 


789 


533 


227 


428 


2 





.286 


614 


046 


845 


766 


0. 


.238 


145 


631 


122 


015 





238 


145 


631 


122 


015 


3 





.200 


705 


485 


068 


072 


0. 


.238 


457 


646 


856 


530 





.238 


437 


505 


168 


542 


4 





.269 


409 


131 


416 


317 


0. 


.238 


433 


043 


073 


649 





.238 


434 


595 


265 


258 


5 





.212 


175 


660 


263 


795 


0. 


.238 


434 


334 


888 


202 





.238 


434 


314 


202 


617 


6 





.261 


216 


282 


820 


282 


0. 


.238 


434 


307 


597 


861 





238 


434 


305 


531 


575 


7 





.218 


319 


995 


918 


563 


0. 


.238 


434 


299 


034 


546 





238 


434 


301 


380 


046 


8 





.256 


437 


590 


806 


900 


0. 


.238 


434 


297 


734 


183 





238 


434 


298 


863 


521 


9 





.222 


142 


795 


720 


768 


0. 


.238 


434 


298 


685 


404 





238 


434 


298 


539 


687 


10 


0. 


.253 


309 


971 


367 


287 


0. 


.238 


434 


298 


885 


696 





.238 


434 


298 


711 


196 


11 





.224 


749 


013 


245 


827 


0. 


.238 


434 


298 


760 


156 





238 


434 


298 


769 


407 


12 





.251 


104 


853 


004 


796 


0. 


.238 


434 


298 


751 


345 





238 


434 


298 


766 


552 


13 





.226 


638 


970 


976 


291 


0. 


.238 


434 


298 


765 


877 





238 


434 




(DO 


311 


14 





.249 


467 


024 


022 


419 


0. 


.238 


434 


298 


763 


853 





238 


434 


298 


763 


230 


15 





.228 


071 


952 


328 


669 


0. 


.238 


434 


298 


762 


980 





238 


434 


298 


763 


331 


16 





.248 


202 


731 


012 


486 


0. 


.238 


434 


298 


763 


382 





238 


434 


298 


763 


332 


17 





.229 


195 


680 


986 


673 


0. 


.238 


434 


298 


763 


343 





238 


434 


298 


763 


330 


18 





.247 


197 


365 


060 


190 


0. 


.238 


434 


298 


763 


322 





.238 


434 


298 


763 


330 


19 





.230 


100 


443 


575 


937 


0. 


.238 


434 


298 


763 


332 





238 


434 


298 


763 


330 


20 


0. 


.246 


378 


830 


994 


286 


0. 


238 


434 


298 


763 


330 


0. 


.238 


434 


298 


763 


330 



exact 0.238 434 298 763 330 0.238 434 298 763 330 0.238 434 298 763 330 



TABLE X. Evaluation of 10" 1 3 F 2 (1, 3, 7; 5/2, 14; 0.99999) with the CNCT. 



n 






t 


3„ 








di 0) ( 


'l,S 


o) 






5i 0) ( 


l,So) 







0. 


.354 


205 


299 


194 


014 


0. 


.354 


205 


299 


194 


014 


0. 


.354 


205 


299 


194 


014 


1 





.227 


102 


649 


597 


007 


0. 


.268 


438 


401 


594 


236 





.268 


438 


401 


594 


236 


2 





.288 


360 


357 


978 


268 


0. 


266 


943 


489 


834 


494 





.266 


943 


489 


834 


494 


3 





.254 


808 


733 


179 


764 


0. 


.267 


121 


290 


068 


426 





,267 


112 


224 


310 


036 


4 





.274 


632 


798 


536 


574 


0. 


.267 


100 


105 


730 


057 





,267 


101 


775 


442 


210 


5 





.262 


289 


292 


919 


201 


0. 


.267 


103 


299 


447 


369 





,267 


102 


948 


107 


378 


6 





.270 


285 


887 


617 


646 


0. 


.267 


102 


742 


191 


670 





,267 


102 


815 


189 


177 


7 





.264 


938 


303 


793 


252 


0. 


.267 


102 


836 


668 


825 


0. 


,267 


102 


824 


285 


564 


8 





.268 


610 


033 


794 


438 


0. 


.267 


102 


822 


243 


079 





,267 


102 


823 


985 


155 


9 





.266 


031 


550 


394 


689 


0. 


.267 


102 


824 


198 


445 





,267 


102 


823 


985 


352 


10 





.267 


878 


112 


865 


182 


0. 


.267 


102 


823 


960 


840 





,267 


102 


823 


984 


751 


11 





.266 


532 


664 


737 


038 


0. 


.267 


102 


823 


987 


269 





,267 


102 


823 


984 


758 


12 





.267 


528 


211 


474 


989 


0. 


.267 


102 


823 


984 


510 





,267 


102 


823 


984 


761 


13 


0. 


.266 


781 


286 


870 


185 


0. 


.267 


102 


823 


984 


786 





,267 


102 


823 


984 


762 


14 


0. 


.267 


348 


762 


560 


130 


0. 


.267 


102 


823 


984 


759 





,267 


102 


823 


984 


762 


15 


0. 


.266 


912 


657 


747 


879 


0. 


.267 


102 


823 


984 


762 





,267 


102 


823 


984 


762 


16 





.267 


251 


339 


082 


039 


0. 


.267 


102 


823 


984 


762 





,267 


102 


823 


984 


761 


17 





.266 


985 


765 


115 


182 


0. 


.267 


102 


823 


984 


761 





,267 


102 


823 


984 


762 


18 





.267 


195 


878 


973 


722 


0. 


.267 


102 


823 


984 


762 





,267 


102 


823 


984 


762 


19 





.267 


028 


262 


510 


152 


0. 


.267 


102 


823 


984 


762 





,267 


102 


823 


984 


762 


20 


0. 


.267 


163 


009 


854 


165 


0. 


.267 


102 


823 


984 


762 





,267 


102 


823 


984 


762 



exact 0.267 102 823 984 762 0.267 102 823 984 762 0.267 102 823 984 762 
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TABLE XI. Evaluation of 1CT 1 3 F 2 (1, 3, 7; 5/2, 14; 1) with the CNCT. 



n 






t 


3„ 










di 0) ( 


'l,S 


o) 






5i 0) ( 


l,So) 







0. 


.354 


212 


896 


979 


703 


0. 


.354 


212 


896 


979 


703 


0. 


,354 


212 


896 


979 


703 


1 





.227 


106 


448 


489 


852 


0. 


.268 


443 


680 


394 


043 





,268 


443 


680 


394 


043 


2 





.288 


366 


524 


620 


961 


0. 


266 


948 


705 


538 


902 





,266 


948 


705 


538 


902 


3 





.254 


813 


300 


376 


035 


0. 


.267 


126 


514 


679 


686 





,267 


117 


448 


402 


341 


4 





.274 


638 


493 


612 


452 


0. 


.267 


105 


329 


111 


381 





,267 


106 


998 


932 


606 


5 





.262 


294 


169 


832 


612 


0. 


.267 


108 


523 


033 


192 





,267 


108 


171 


668 


587 


6 





.270 


291 


370 


710 


880 


0. 


.267 


107 


965 


739 


446 





,267 


108 


038 


742 


291 


7 





.264 


943 


330 


016 


988 


0. 


.267 


108 


060 


223 


478 





,267 


108 


047 


839 


255 


8 





.268 


615 


409 


392 


370 


0. 


.267 


108 


045 


796 


597 





,267 


108 


047 


538 


821 


9 





.266 


036 


655 


402 


122 


0. 


.267 


108 


047 


752 


131 





,267 


108 


047 


539 


018 


10 


0. 


.267 


883 


429 


838 


931 


0. 


.267 


108 


047 


514 


503 





,267 


108 


047 


538 


417 


11 





.266 


537 


813 


952 


027 


0. 


.267 


108 


047 


540 


935 





,267 


108 


047 


538 


424 


12 





.267 


533 


494 


713 


522 


0. 


.267 


108 


047 


538 


176 





,267 


108 


047 


538 


427 


13 





.266 


786 


462 


107 


999 


0. 


.267 


108 


047 


538 


452 





20 / 


1 HQ 


U4/ 


£ Q Q 


428 


14 





.267 


354 


025 


525 


918 


0. 


.267 


108 


047 


538 


425 





,267 


108 


047 


538 


428 


15 





.266 


917 


848 


923 


282 


0. 


.267 


108 


047 


538 


428 





,267 


108 


047 


538 


428 


16 





.267 


256 


589 


412 


242 


0. 


.267 


108 


047 


538 


428 





,267 


108 


047 


538 


428 


17 





.266 


990 


966 


387 


224 


0. 


.267 


108 


047 


538 


428 





,267 


108 


047 


538 


428 


18 





.267 


201 


121 


176 


708 


0. 


.267 


108 


047 


538 


428 





,267 


108 


047 


538 


428 


19 





.267 


033 


470 


369 


207 


0. 


.267 


108 


047 


538 


428 





,267 


108 


047 


538 


428 


20 


0. 


.267 


168 


246 


683 


931 


0. 


267 


108 


047 


538 


428 


0. 


,267 


108 


047 


538 


428 



exact 0.267 108 047 538 428 0.267 108 047 538 428 0.267 108 047 538 428 



TABLE XII. Evaluation of 10~ 5 Ei= ( 2 ' + 0.9999 x 0.7) /i ; (1) (i0.7) with the CNCT. 



n 






S 


n 








a 


'n" 1 (l, So 


) 






8 


< o) 


, S u 


) 







-0. 


206 


084 


520 


894 


668 


-0. 


.206 


084 


520 


894 


668 


-0. 


206 


084 


520 


894 


668 


1 


-0. 
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